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ABSTRACT 

We present ultraviolet, optical and near-infrared observations of the interacting tran- 
sient SN 2009ip, covering the period from the start of the outburst in October 2012 
until the end of the 2012 observing season. The transient reached a peak magnitude 
of My=-YJ .1 mag, and with a total integrated luminosity of 1.9xl0 49 erg over the 
period of August-December 2012. The lightcurve fades rapidly, dropping by 4.5 mag 
from the V band peak in 100 days. The optical and near infrared spectra are domi- 
nated by narrow emission lines with broad electron scattering wings, signaling a dense 
circumstellar environment, together with multiple components of broad emission and 
absorption in H and He at velocities between 0.5-1.2 xlO 4 krns" 1 . We see no evidence 
for nucleosynthesized material in SN 2009ip, even in late-time pseudo-nebular spectra. 
We set a limit of <0.02 M© on the mass of any possible synthesized 56 Ni from the late 
time lightcurve. A simple model for the narrow Balmer lines is presented, and used to 
derive number densities for the circumstellar medium of between ~ 10 9 — 10 10 cm~ 3 . 
Our near-infrared data does not show any excess at longer wavelengths, and we see 
no other signs of dust formation. Our last data, taken in December 2012, shows that 
SN 2009ip has spectroscopically evolved to something quite similar to its appearance 
in late 2009, albeit with higher velocities. It is possible that neither of the eruptive 
and high luminosity events of SN2009ip were induced by a core-collapse. We show 
that the peak and total integrated luminosity can be due to the efficient conversion of 
kinetic energy from colliding ejecta, and that around 0.05-0.1 M of material moving 
at 0.5-1 xlO 4 kms _1 could comfortably produce the observed luminosity. The ejection 
of multiple shells, lack of evidence for nucleosynthesied elements and broad nebular 
lines, are all consistent with the pulsational-pair instability scenario. In this case the 
progenitor star may still exist, and will be observed after the current outburst fades. 
The long term monitoring of SN 2009ip, due its proximity, has given the most exten- 
sive data set yet gathered of a high luminosity interacting transient and its progenitor. 
It is possible that purported Type Iln SNe are in fact analogues of the 2012b event 
and that pre-explosion outbursts have gone undetected. 

Key words: stars: massive — supernovae: general — supernovae: individual 
(SN2009ip) 



1 INTRODUCTION 

Stars with a zero-age main sequence (ZAMS) mass of greater 
than ~8 Mq end their lives as core-collapse supernova 
(CCSNe). The differing masses, metallicities and structures 
which these massive stars can have at the point of collapse 
give rise to a tremendous diversity of SN types. Among the 
most peculiar a nd extreme are the Type Iln (narrow) SNe 
|Schlegel|[T990l '). which are typified by narrow (-100-500 
kms^ 1 ) emission lines of H, which arise from the collision of 
the SN eje cta with pre-existing dense ci rcumstellar material 
(CSM; see IChevalier fc Franssonl (|l994l ) for an overview of 
the theory of ejecta-CSM interaction). 

Type Iln SNe are intrinsically rare, with relative rates 
(i.e. rates as a fraction of the total nu mber of CCSNe) es- 
timated to lie betw een 4 to 9 per cent l|Smartt et al.l [20091 : 
ISmith et all l2011bl ). Their rarity is compounded by their 
diverse nature; with well studied e xamples including long - 
lived events (SNe 2005ip and 2006id; IStritzinger et alfeoia) 
events which show Wolf -Rayet features (SN 1998S; Liu et all 
2000 , iFassia et al.ll2001f ) , and even a related family of events 
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which show interaction with a H e-rich rather than H-rich 
CSM l|Pastorello et al.ll200ll200sl ). 

In some cases, Type Iln SNe have been linked to Lu- 
minous Blue Variables (LBVs). LBVs are characterized by 
having bright absolute magnitudes (Mtoi < —9.5), irregular 
eruptions where they brighten by more than one magnitude, 
and a s urface temperature during qui escence of between 12- 
30 kK ^Humphreys fc Davidsonll 19941 ). As Type Iln SNe re- 
quire a dense circumstellar environment prior to explosion, a 
connection to LBVs - which shed large amounts of material 
into their CSM - seems natural. 

The first connection between LBVs and Type Iln SNe 
was that eruptions of the fo rmer are often confu sed with 
the latter - a fact which led IVan Dvk et al. I (|2000l ) to coin 
the term "supernova impostors" to describe extreme out- 
bursts of massive stars. While SN impostors share the key 
characteristics of Type Iln SNe, viz. narrow emission lines 
signaling interaction, the crucial difference is that a SN 
impostor is a non-terminal event, whereas the progenitor 
star will not survive a SN explosion. The boundary be- 
tween Type Iln SNe and impostors has been m uch de- 
bated in the literature (see, f or example, SN 1961V ijZwickv 
I l96l|) , whi ch was revisited bvlKochanek. Szczygiel fc Stanek 



and IVan Dvk fc Mathesonl 2012 



SN 1994W-like objects. ISollerman et al 



120091; iKankare et alJliuli i 



and the debated 
ll998l ; lDessart et all 
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Direct evidence for a n LBV- Type Iln SN link was found 
bv lGal-Yam et all (|2007l ). who identified a source coincident 
with SN 2005gl in pre-explosion imaging, which was subse- 
quently seen to have disap peared after the SN had faded 
l|Gal-Yam fc LeonardlixH ) . The bright absolute magnitude 
of the progenitor for SN 2005gl could only be matched by 
an LBV. Similar resu lts from archival im aging have been 
found for SN 2010jl bv lSmith et~ail (|2011al ). A physical con- 
nection betwe e n Typ e Iln SNe and LBVs was also made by 
iTrundle et alj l|200a ). who observed multiple narrow absorp- 
tion components for H and He in the spectrum of the Type 
Iln SN 2005gj, which they attributed to var iable mass loss 
in an LBV progenitor (although we note that lAldering et al.l 
(2006) have proposed an alternative scenario, where SN 
2005gj was a Type la SN exploding in a dense CSM). 

In August 2009, SN 2009i p was discove red close to NGC 
7259 by the CHASE project (|Maza et al.ll2009h . While the 
object was given an SN designation, it was rapidly real- 
ized th^Jt_was no^a genuine SN, but rather an impos- 
tor l|Smith et al.ll2010al ; iFolev et al.ll201ll ). Fortuitously, the 
site of SN 2009ip had been observed with the Hubble Space 
Tel escope prior to the 2009 outburst, enabling ISmith et al.l 
and lFolev et alj to identity a coincident point source with an 
absolute magnitude of F606W=-9.8; and consistent with a 
massive star of ~50-60 Mq . While the August 2009 eruption 
faded, SN 2009ip continu ed to suffer irregu l ar ou tbursts over 
the following three years. IPastorello et alj |2013l . henceforth 
P13) presented coverage of these eruptions, with a unique set 
of spectroscopic and photometric data, making SN 2009ip 
one of the best intensively monitored extragalactic tran- 
sients to date. Remarkably, P13 found evidence for material 
accelerated to high velocities (~ 10 4 kms^ 1 ) in September 
2011. The detection of fast material over one year prior to 
the 2012a and 2012b events - and at a point when SN 2009ip 
had not undergone core-collapse - shows that high velocities 
cannot be taken as proof of a CCSN. 

In July-August 2012, SN 2 009ip started a new phase 
of activity, as first reported by iDrake et alj (|2012h - After 
an initial outburst which reached a peak absolute mag- 
nitude of V=-14 and lasted until around September 25 
(termed the "2012a" event), SN 2009ip rapidly rebright- 
ened to reach an unprecedented pea k magnitude of J ? =-18.5 
ma g (the "2012b" out burst). P13, iMauerhan et al.l \20l3 ) 
and lPrieto et al.l (|2013l ') all published early photometry and 
spectroscopy of the eruption, which showed both signs of 
interaction (narrow lines) and broad absorptions in H and 
He, with velocities of up to 15000 kms -1 . 

In this paper we present high cadence ultraviolet (UV), 
optical and near-infrared (NIR) photometry of SN 2009ip, 
together with optical and NIR spectroscopy, covering the 
period from the peak of the 2012b event, until the end of 
the observing season in December 2012. The rest of the pa- 
per is organized as follows: in Section [2] we describe the 
observations obtained, and the reduction techniques used. 
In Sections [3] and [4] we characterize the photometric and 
spectroscopic evolution of the SN, and derive a bolometric 
lightcurve. Section \S\ is devoted to simple modelling of the 
narrow H emission seen in SN 2009ip, which we use to de- 
rive densities for the CSM. Finally, in Section |6j we attempt 
to use our observational data and modelling to understand 
the evolution of SN 2009ip, and discuss whether the out- 
bursts seen in late 2012 represent the final death of the star 



as a SN, or rather an example of a non-terminal outburst. 
In all of the following, we adopt the foreground reddening 
and distance towards SN 2009ip quoted by P13 and others; 
^i=31.55 mag, A.r=0.05 mag. The distance adopted i s also 
in agreement with that found bv lPotashov et al.l l|2012l ) from 
applying a geometric technique (the "Dense Shell Method" ) 
to SN 2009ip itself. 



2 OBSERVATIONS AND DATA ANALYSIS 

Optical spectroscopic follow-up of SN 2009ip was chiefly 
obtained with the New Technology Telescope (NTT) + 
EFOSC2, as part of the Public ESO Spectroscopic Sur- 
vey of Transient Objects (PESSTOfl The PESSTO data 
were supplemented with data from the Telescopio Nazionale 
Galileo (TNG) + DOLORES, and the Australian National 
University 2.3m telescope + WiFeS. Details of the instru- 
mental configuration used for each spectrum are listed in Ta- 
bleQ] All PESSTO data are also immediately publicly avail- 
able from the ESO archive; high level data products will be 
available in June- July 2013. All spectroscopic data presented 
in this paper will be available for download from the Weiz- 
mann Interactive Super nova data REPository (WISEREP; 
lYaron fc Gal-Yaml l2012'PI upon acceptance. 

All long slit spectra were reduced in the standard fash- 
ion; bias and overscan subtracted, flat-fielded using normal- 
ized lamp flats, wavelength-calibrated to an arc lamp, and 
flux-calibrated to a spectophotometric standard star. In the 
case of the EFOSC2 spectra obtained by PESSTO, these 
steps were carried out within the PESSTO pipelin^fl, while 
for the DOLORES spectra the reduction was performed 
within iraiQ The PESSTO data were corrected for telluric 
absorption by subtracting a model spectrum of the telluric 
bands, while the DOLORES spectra were corrected using a 
smooth spectral standard. For spectra where there were a 
signific ant number of cosm ic ray hits, the LACOSMIC algo- 
rithm l|Van DokkumllioOlh was used to identify and mask 
cosmic rays in the two-dimensional frames before extrac- 
tion. The resolutions of the slit spectra were checked from 
the FWHM of narrow night sky lines, and were found to 
be ~15 A and ~10 A for the EFOSC2 and DOLORES data 
respectively. These values are in good agreement with the ex- 
pected resolution from the grating characteristics (reported 
in Table [TJ and correspond to a resolution in velocity space 
(at the wavelength of Ha) of ~500-700 kms -1 . 

Some optical spectra of SN 2009i p were obtained with 
the Wide Field Spectrograph (WiFeS; iDopita I [20071 . l2010h 
on the Australian National University (ANU) 2.3 m tele- 
scope at Siding Spring Observatory in northern New South 
Wales, Australia. Observations were obtained with either 



1 www.pessto.org 

2 http:/ /www. weizmann.ac.il/astrophysics/wiserep/ 

3 The PESSTO pipeline has been developed by S. Valenti, and 
comprises a set of PYTHON scripts which call PYRAF tasks to reduce 
EFOSC2 and SOFI data. 

4 I RAF is distributed by the National Optical Astronomy Observa- 
tory, which is operated by the Association of Universities for Re- 
search in Astronomy (AURA) under cooperative agreement with 
the National Science Foundation. 
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the B3000/R3000 gratings, which provide wavelength cov- 
erage from 3500 to 9600 A at a resolution of ~2 A, or the 
U/B/R/I7000 gratings which cover 3500 to 9100 A with ~1 
A resolution. As WiFeS is an integral-field spectrograph, the 
data were not reduced using the same techniques as for the 
rest of the optical spectra. Instead, data cubes for WiFeS ob- 
servations were produced using the pywifes softwar«[il and 
final spectra were obtained with a PSF-weighted extraction 
routine. 

NIR spectroscopy was obtained with the NTT+SOFI 
and LBT+LUCIFER as detailed in Tabled Solar analogues 
at a similar airmass were observed either before or after SN 
2009ip, to facilitate the removal of the strong telluric ab- 
sorptions which lie between 1 and 2u.m. For two epochs 
of NIR spectra (on Nov 13 and Dec 4), a flux standard 
was observed on the same night and used to calibrate the 
spectra, for all other data the telluric standard was also 
used to flux calibrate SN 2009ip. The LUCIFER spectrum 
was reduced with the ad-hoc pipeline of the LBT Spec- 
troscopic Reduction Center of IASF-Milano, which includes 
bad pixel, dark and flat-field correction, optimal extraction, 
sky-subtraction, and wavelength and flux calibration. 

Optical photometry was obtained with the NTT + 
EFOSC2 + t/#640,B#639y#641,i?#642/#705 (on the 
same nights as the EFOSC2 optical spectroscopy); with the 
Panchromatic Robotic O ptical Monitoring a nd Polarime- 
try Telescope (PROMPT; iReichart et alj|2005h telescopes + 
BVRI on Cerro Tolodo, and with the Liverpool Telescope 
+ RATCAM + u'BVr'i'. In addition, regular photometry 
of SN 2009ip was obtained with the Ultra Violet and Optical 
Telescope (UVOT) onboard the Swift satellite (Programme 
IDs 31486, 32579; P.I. Margutti). The PROMPT and RAT- 
CAM images were reduced automatically (trimmed, bias- 
subtracted and flat-fielded) by their respective pipelines. 
Significant fringing was evident in the LT RATCAM i' im- 
ages, and so a fringe frame was constructed from deep expo- 
sures of blank fields, scaled to match the science images, and 
subtracted from them. The EFOSC2 images were reduced 
using the PESSTO pipeline, again subtracting a smoothed 
fringe frame in I. 

For the PROMPT and RATCAM data, point-spread 
function (PSF) fitting photometry was performed using the 
SNOOPY package within IRAF. SNOOPY is a PSF-fitting pack- 
age using the routines in the daophot package, but specifi- 
cally designed for SN photometry, and including several ad- 
ditional functions such as artificial star tests to estimate 
photometric errors. The zero point for each image was de- 
termined from aperture photometry of the sequence stars 
listed in P13, while the same colour terms as used by P13 
were applied to the instrumental magnitudes. The smaller 
field of view of EFOSC2, coupled with the short exposure 
times necessary to avoid saturating the bright target, re- 
sulted in only a single reference star bring present in the 
images obtained in October. As this is insufficient to build a 
PSF for each frame, we performed aperture photometry on 
these data. As for the PROMPT and LT data, the zeropoint 
for each night and filter was determined from local sequence 
stars, while colour terms for EFOSC2 were evaluated from 
observations of Landolt standard fields on three photomet- 



http: / /www.mso . anu.edu.au/pywifes/ 



ric nights. All optical photometry, calibrated to the Landolt 
photometic system is reported in Tabic [Bll 

The Swift UVOT frames were reduced using tools 
within HEAsofifl Individual images for each epoch were first 
co-added, before apert ure magnitudes we re measured follow- 
ing the prescription of lPoole et al.l (|2008T ). To transform the 
UVOT ubv magnitudes from the instrumental system into 
standard UBV Johnson magnitudes, shifts of AU = 0.270, 
AB = 0.018 and AV = -0.042 mag were applied. These 
shifts were computed from the magnitudes of the reference 
stars in the SN field, and are consistent with those reported 
in P13. These corrections have been applied to all UVOT 
magnitudes and the resulting magnitude values are reported 
in Table lBTI The uvw2, uvm2, uvwl magnitudes are reported 
in Table IB2I in the UVOT photometric system. 

NIR imaging of SN 2009ip in JHKs were obtained with 
the NTT + SOFI, the LBT + LUCIFER, and with the 
Nordic Optical Telescope + NOTCAM. Multiple, dithered, 
on-source images were taken in each filter, these images were 
then flat fielded and median combined to create a sky frame. 
The sky frame was then subtracted from each of the individ- 
ual images, which were then aligned and coadded. For the 
LUCIFER imaging, these steps were carried out with the 
pipeline at the LBT Imaging Data Centre at the Observa- 
tory of Rome. Aperture photometry was performed on SN 
2009ip, taking care to avoid any flux from the nearby fore- 
ground star to the north. Zeropoints for each image were 
determined from aperture photometry of the same sequence 
stars as used for the optical photometry, calibrating them 
to their 2MASS JHK magnitudes. While the K and Ks fil- 
ters are slightly different, this has a negligible effect on our 
photometry. NIR magnitudes for SN 2009ip are reported in 
Table |H 



3 SPECTROSCOPIC DATA 

The full set of spectra collected for SN 2009ip is presented 
in Tables [T] and [2] In general, the spectra are blue at early 
times, and are dominated by narrow emission lines from H 
and He I. As the SN evolves, broad (~ 10 4 kms _1 ) absorp- 
tions appear in the Balmer lines, while the [Ca n] NIR triplet 
emerges in emission, along with Na I D. We do not observe 
a transition to a nebular SN spectrum; indeed the last spec- 
trum in our data set appears to resemble most closely the 
spectrum of SN 2009ip from September 2009, with the same 
lines present at both epochs, albeit at higher velocities in 
2012. In Sections 13. II and 13. 21 we describe in detail the spec- 
troscopic evolution of SN 2009ip. 



3.1 Optical spectoscopic evolution 
• Sept 22 - Sept 23 

The first spectrum in our sequence (from Sept 22) 
as shown in Fig. [1] is contemporaneous with the 
TNG+DOLORES spectrum presented in P13; coinciding 
with the end of the 2012a eruption, and immediately prior 



5 available from the NASA High Energy Astrophysics Science 
Archive Research Center 
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3000 4000 5000 6000 7000 8000 9000 10000 

Rest wavelength (A) 

Figure 1. The sequence of optical spectra of SN 2009ip, as detailed in Table [T] Note that the high resolution spectrum obtained on 23 
Oct is not included here, but is shown separately in Fig. [3] The wavelength of residuals from the telluric correction are indicated with ®. 
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Table 1. Optical spectroscopy of SN 2009ip, together with the 
nominal resolution for each grating. 
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Table 2. NIR spectroscopy of SN 2009ip 
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Figure 2. The region of Ha in our Oct 8 spectrum of SN 2009ip, 
showing the shallow broad absorption. In red is a blackbody fit 
to the continuum, while in green is the same blackbody, with a 
Lorentzian line of FWHM 17 A and a Gaussian absorption with 
FWHM 270 A superimposed. The Lorentzian plus continuum is 
an excellent fit to the red wing of Ho (with the exception of He 
A6678), whereas in the blue wing there is an apparent absorption 
which is fitted by the Gaussian. As discussed in the text, it ap- 
pears that the absorption may have two components, although a 
single Gaussian (as shown here) also gives a reasonable fit. The 
second x axis shows the velocity with respect to the rest wave- 
length of Hce. 
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to the start of the 2012b event. As this spectrum is effec- 
tively identical to one obtained on the following night (Sept 
23), these two spectra will be discussed together. The most 
prominent features in the spectra are the Balmer lines. At 
the rest wavelength of Ha, SN 2009ip has a narrow (FWHM 
~500 kms -1 ) emission component superimposed on broad 
emission (FWHM ~6000 kms -1 ). A P-Cygni absorption 
trough is also present with a minimum at -7500 kms -1 with 
respect to the rest wavelength of Ha. This absorption is 
weaker than the broad emission. The higher Balmer lines, 
notably H/3 and H7, have similar multi-component profiles, 
but with a stronger broad absorption compared to the broad 
emission. The minimum of the broad absorption in H/3 is also 
at a lower velocity than that seen in Ha. A "shoulder" of 
emission is also visible in the blue wing of H/3 and H7 (at 
a velocity of -1800 kms -1 ), which is less visible in Ha. The 
He 1 5876 A line which appears later is absent at this phase. 
Metals are present as narrow emission lines, including the 
blended Na I D lines at 5890 A and 5895 A, weak Ca 11 NIR 
triplet lines and numerous Fe 11 lines such as Multiplets 27 
and 42. 

• Oct 8 - Oct 15 

The Oct 8 and Oct 9 spectra were obtained when SN 
2009ip was at the peak of the 2012b event. The spectra 



have a blue continuum, and are dominated by narrow (un- 
resolved, <700 kms -1 ) Balmer emission. The Paschen se- 
ries are also visible at wavelengths >8500 A although with 
much weaker emission. The narrow central component of 
Ha, together with the broad base of the line are well fit by a 
Lorentzian with a FWHM 750km s -1 . The broad absorption 
is still present, but much weaker than in the previous epoch. 
We show in detail the Ha line profile in Fig. [2] together with 
a Lorentzian fit, clearly showing the broad shallow absorp- 
tion in the blue. Tantalizingly, the profile appears to resem- 
ble the multi-component absorption we see at later epochs; 
and we obtain a best fit to the absorption component with 
two Gaussians with minima at -6000 and -12000 kms -1 . 
However, as the feature is weak we cannot formally exclude 
a single broad absorption. He I A5876 is the strongest line 
after the Balmer series, it has a similar Lorentzian profile to 
Ha (narrow emission on a broad base, together with a hint 
of broad absorption in the blue). Other weaker He I lines are 
detected at 4471, 4713, 5015, 6678, 7065 A and tentatively 
at 7276 A. We see no lines from He n (for example at 4686 
A), nor of the Ca NIR triplet or Na I D. 

The next spectra in our sequence are from Oct 14 and Oct 
15. These spectra are similar to those obtained on Oct 8 and 
9, but with a double absorption in the Balmer lines. Fitting 
two Gaussians to the absorptions in Ha, we find velocities 
of -8000 and -12500 kms" 1 . Interestingly, while H/3 and H7 
appear to have the same double absoption as Ha, the lower 
velocity component is at only ~4000 kms -1 . The higher ve- 
locity component is difficult to measure accurately due to 
the steep continuum, but it too seems to be lower than in 
Ha. He I 5876 A has a similar narrow Lorentzian emission 
component to the preceding epoch, but now also has a clear 
broad absorption centered on -12500 kms -1 . This absorp- 
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tion does not appear to have the two-component structure 
seen in the Balmer lines, but is instead well fitted by a sin- 
gle Gaussian absorption. The broad absorption is not seen 
in the other He I lines, but as these are much weaker than 
He I 5876 A this is perhaps not surprising. While an alterna- 
tive source for the broad absorption could be Na, this would 
imply a velocity of 13500 kms -1 , which is uncomfortably 
higher than Ha. 

• Oct 20 - Oct 23 

On Oct 22 we obtained spectra of SN 2009ip from three 
different telescopes, all of which are similar to that obtained 
two days earlier on Oct 20. At this phase, the narrow emis- 
sion on a broad base in Ha can no longer be fitted with 
a single Lorentzian profile. Instead, we find that we need 
two components to fit the emission; one narrow (FWHM 
~800 kms" 1 ) and one broad (FWHM ~10000 kms" 1 ). The 
broad emission is not symmetric, as a "shoulder" of emis- 
sion is present in the blue wing of Ha which is not present 
in the red. While the broad absorption can be reasonably fit 
with a single Gaussian at -10000 kms -1 , if we take the po- 
sition of the two absorption components from the previous 
epoch, then we can see two notches still present at veloc- 
ities of ~8000 and 12000 kms" 1 . H/3 has a similar broad 
absorption to Ha, but with a much less pronounced broad 
emission, while in H5 the broad emission is absent. For He I, 
the narrow emission is much weaker than in previous epochs, 
and there is no strong broad emission, while we still see a 
strong absorption at -11000 kms -1 in the 5876 A line. At 
this phase, we also see the emergence of a broad emission 
feature at ~4500 A with a FWHM of ~10000 kms -1 . It is 
difficult to judge the strength of the 4500 A feature, as it 
is bracketed by the broad absorption of H/3 and H<5, which 
serve to lower the apparent level of the continuum. While 
this feature was present during Oct 14 and 15, it was much 
weaker, and was plausibly associated with the He 4471 A 
line. However, for the 4500 A line to be caused by He I 
on Oct 22 would require the peak of the emission to be red- 
shifted by ~2000 kms" 1 . More problematically, it is difficult 
to see how we could get emission in the 4471 A line, without 
corresponding emission in the 5876 A He line. While there 
are no strong, low-ionization state lines arising from any in- 
termediate mass elements around 4500 A we note that in 
Type II SNe an emission feature at this wavelength is of- 
ten seen at early times, arising from multiple blended Fe I 
and Fe II lines such as multiplets 37 and 38. The Fe which 
gives rise to this in SNe is primordial, so it is entirely possible 
that we could see a similar feature in a massive non-terminal 
eruption of an LBV. We also see two emission lines (FWHM 
~50 A) in the extreme blue end of the spectrum, at 3490 and 
3650 A. The 3650 A feature may be caused from the com- 
bined flux of the higher order Balmer lines (10-2 transitions 
and higher), as the series converges to this wavelength, and 
the blue edge of the line corresponds to the wavelength at 
which the separation between the lines becomes lower than 
the instrumental resolution of EFOSC2 (for comparison we 
can clearly resolve the Hj=g_2 line in the same spectrum). 
The identity of the 3490 A line is uncertain, although as for 
the ~4500 A bump, it may be attributable to blended Fe 
emission. 

We obtained a higher resolution spectrum of SN 2009ip on 
October 23 with WiFeS using the R=7000 gratings (FigEJ). 
We measured a FWHM of the narrow Lorentzian emission in 



Ha of 266±13 kms" 1 , 243 kms" 1 in H/3 and 242 kms" 1 in 
H7. As the FWHM is larger than the instrumental resolution 
of WIFES (45km s" 1 at Ha), we adopt their mean (250±11 
kms" 1 ) as the velocity of the CSM. The 3650 A feature is 
indeed seen to be comprised of a large number of Balmer 
lines, as the blue edge of the bump as seen in the spectra 
from Oct 22 are resolved out into a blend of the Balmer 
lines with upper levels between J=13 and J=22. The other 
striking feature is the appearance of a large number of weak, 
narrow emission lines from Fe. These features are not seen 
in the lower resolution data from the preceding night, but 
were present in spectra of earlier eruptions obtained with 
XShooter by P13. Unfortunately we do not observe the au- 
roal O line at 5007 A precluding a direct estimate of the 
metallicity of the environment. 
• Oct 30 - Nov 7 

In the spectra obtained on Oct 30 and Nov 04, Ha has 
quite distinct narrow and broad emission features. The for- 
mer has a FWHM of ~700 kms" 1 , while the latter has a 
FWHM ~10000 kms" 1 . The broad absorptions seen at ear- 
lier epochs are still present, with multiple components at 
12500 and 10000 kms -1 . We note however that the broad 
absorptions are weaker than the corresponding broad emis- 
sion by a factor of a few. The blue shoulder in Ha which 
was beginning to appear in the Oct 22 spectrum is more ev- 
ident. If caused by an additional emission component, then 
it consists of material moving at a velocity of -3000 kms 
with respect to the narrow component, and with a FWHM 
of ~1500 kms" 1 . Alternatively, the feature could be viewed 
as an absorption centered on -1500km s , although this is 
perhaps less likely, as it would imply that the strong broad 
emission is intrinsically quite asymmetric. H/3 and H7 do 
not show as strong narrow emission when compared to Ha. 
The broad features in H/3 and H7 are also quite peculiar. For 
both lines, the broad emission seen strongly in Ha is absent, 
while the broad absorption is much more pronounced. We 
see multiple absorptions, with the highest velocity minimum 
at -12500 kms^'in H/3, but with the bulk of the absorp- 
tion at -8000 kms" 1 . We also see for the first time in the 
2012b event clear absorptions from ejecta at quite low ve- 
locities, with the first minimum in H/3 at only -2000 kms -1 . 
The strong absorption at 5660 A which we attribute to He I 
with a velocity of -llOOOkms" 1 is still present, without any 
clear corresponding broad emission. However, the narrow 
He 1 5876 A line has now vanished, and has been replaced 
with a re-emergent Na I D line at 5890,5895 A suggesting 
that from this point onwards we cannot be confident that 
the broad feature is solely caused by He. We also observe 
the re-appearance of narrow Fe lines in the spectrum of SN 
2009ip, such as the 5167 A line from Multiplet 42. 

The spectrum of SN 2009ip from Nov 7 is similar to those 
discussed in the preceding paragraphs. Ha now has a clear, 
three component absorption, with minima at -12500, -10000 
and -6500 kms" 1 . The absorption in H/3 now begins even 
more sharply at a velocity of only -1000 kms -1 . Most in- 
terestingly however, is the confirmation of a two-component 
absorption in the line at 5660 A. If the line is still dominated 
by He, then the velocities of the two absorptions are 6000 
and llOOOkms -1 , while if Na is dominant, the velocities are 
slightly higher at 12000 and 7000 kms -1 . Only the Na I D 
narrow emission is present (albeit weakly), but it is proba- 
ble that the broad emission is a combination of both He and 
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Figure 3. The highest resolution spectrum of SN 2009ip in our data set, taken on 23 Oct with the U7000, B7000, R7000 and 17000 
gratings on WiFeS, showing the forrest of narrow Fe emission lines. 
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Figure 4. The evolution of the line profiles for selected elements in SN 2009ip. All panels are on a common velocity scale, centered on 
the rest wavelength of the line in question, with units of 10 3 kms -1 . From top to bottom, the rows correspond to: pre-2012a outburst 
(Sept 22), peak 2012a outburst (Oct 9), first "bump" in the lightcurvc (Nov 7), approximately two months past peak (Dec 3), and the 
latest spectrum (Dec 20). All plots are on a linear flux scale, normalized to the peak of the line in question, with the exception of the 
column on the far right, which shows the Ha line on a logarithmic scale to emphasise the structure at the base of the line. 



Na. At this epoch, we also start to see the emergence of a 
bump around 8500 A which will develop into the Ca NIR 
triplet in later spectra. In the blue part of the spectrum, 
while line identifications are difficult given the complicated 
velocity structure of the outburst, a broad emission line at 
~4100 A appears to be present, which we tentatively iden- 
tify as a blend of Fe n lines. 

• Nov 12 - Nov 20 



By Nov 12, the "shoulder" in the blue wing of Ha has 
begun to fade. We again fit the emission in Ha with a narrow 
and a broad Gaussian, the former is centered on the rest 
wavelength, while the latter is shifted by 1000 kms -1 to 
the red, and has a FWHM of ~10000 kms -1 . The majority 
of the flux in Ha is now emitted by the broad component, 
while the broad absorptions (which are by now extremely 
shallow) persist at -12500, -10000 and -6500 kms" 1 . The 
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overall shape of the absorption in Ha is flat and shallow, a 
shape which is mirrored in H/3. We also see the emergence 
of two emission lines at 5021 and 7075 A. While both these 
features are weak (and so their centers are uncertain), if 
these are the He I 5015 and 7065 A lines, then they are 
both systematically redshifted by several hundred kms -1 . 
This is in contrast to Ha, which is within ~1 A of its rest 
wavelength. We also see what appears to be a broad bump in 
the spectrum at ~5250 A, to the red of the probable He line. 
Superimposed on this bump are numerous weak Fe emission 
lines, indicating that this could be from Fe. In the probable 
He/Na blend at 5660 A, the redder absorption component 
which emerged on Nov 7 is now comparable in strength to 
the other absorption. The Ca NIR triplet is clearly emerging 
in the red part of the spectrum, with high velocities and an 
apparent P-Cygni profile. 

On Nov 20, the complex Ha structure persists with ab- 
sorption up to -12500 kms -1 , H/3 is largely unchanged from 
the preceding spectrum. For H5, the expected absorption is 
masked by the Fe II line at 4100. The narrow He emission 
grows in intensity, as does the Ca triplet. Most interestingly, 
the possible Na/He feature at ~5700 A now has two clearly 
separated absorption minima which are well fitted by Gaus- 
sians, one centered at -13000 kms -1 and one at -6500 krns" 1 . 
The narrow emission is still located at the expected position 
for the unresolved Na I D lines. An apparent absorption fea- 
ture emerges at 5123 A, although at this phase it is very 
difficult to locate the continuum, and we suggest that this is 
only an apparent absorption caused by the strengthening of 
the probable Fe bump between ~5150 - 5550 A, which was 
first seen in the Nov 12 spectrum. 

• Dec 3 - Dec 11 

By December 3, the broad base of the Ha line has become 
quite asymmetric, with significantly more flux in the red 
wing compared to the blue. The narrow emission in Ha and 
H/3 are markedly weaker by this epoch, while H 8 is nearly 
completely gone. Interestingly, the highest velocity absorp- 
tion in Ha, at ~12,000 kms -1 appear to be weakening. Ca 
NIR continues to strengthen, with a clear absorption which 
has a minimum at -3500 kms -1 relative to the bluemost 
component of the the triplet (8498A) . We see the emergence 
of narrow emission features at 6159 and 7313 A which are 
probably Fe. 

On December 11, the spectrum is quite similar to that 
from Dec 3, with the most significant changes being the 
appearance of [Ca n] at 7219,7323 A and the re-emergence 
of strong Paschen lines with a FWHM of ~1500 kms - . 

• Dec 20 

The last optical spectrum we obtained of SN 2009ip was 
on Dec 20. We fit the profile of Ha with a narrow Lorentzian 
emission with FWHM 1200 kms" 1 , and a broad Gaussian 
emission, FWHM 8500 kms -1 . The centre of the broad com- 
ponent is offset by 800 kms - to the red to reproduce the 
asymmetric emission. The absorption component is weaker 
than at previous epochs, and has a blue edge at -11,000 
kms -1 . The [Ca n] emission at 7219,7323 A continues to 
grow stronger, becoming comparable in strength to the He 
I 7065 A line. Meanwhile the P-Cygni profile in the Ca NIR 
triplet has disappeared, to be replaced with a line purely 
in emission, which resembles that of a late-time nebular SN 
line. What is most striking about the final spectrum of SN 
2009ip is the resemblance to the spectra of the transient 



taken three years previously in 2009. A comparison is shown 
in Fig[5] All the lines present in September 2009, are also 
visible in the Dec 2012 spectrum, with the only difference 
being much higher velocities in the latter. The continuum 
is somewhat redder in 2012. H, He, Ca, Na and Fe are all 
present, but we see no evidence for freshly nucleosynthesised 
materials. There is no sign of [O i] at 6300,6364 A, and al- 
though the emission could be speculated to fill-in some of 
the absorption in Ha, we see the same profile in H/3 and Ha, 
indicating that this is not the case. Strong Mg i] 4571 A does 
not appear to be present either, although the presence of Fe 
emission in this region of the spectrum makes this a less 
secure result. However, the lack of clear Mg i] 4571 A is in 
accord with the absence of the 1.5(J. line in the NIR. We also 
see the emergence of a line at ~7710 A, from December 11 
onwards. While the residual from the telluric A-band con- 
taminates this region of the spectrum, the line does appear 
to be real, however, its wavelength is offset by -60 A from O I 
7774. The continuing presence of strong He I emission at this 
epoch is striking, and could be further evidence that inter- 
action is a significant source of energy, as photo-ionization 
will probably be too weak to support these lines. The strong 
He lines may also suggest an overabundance of He, although 
detailed chemical modeling is required to confirm this. 

In Fig.[S]we compare the spectral evolution of SN 2009ip 
to that of a sample of interacting SNe, ranging from the 
weakly interacting Type II SN 2007pk, t hrough to the well 
studied Type Hn SNe 1998S and 201 0jl (|Fassia et all 1200 ll ; 
IZhang et al.ll2012l ; llnserra et alj2012| ). Spectra of SNe 1998S 
and 2010jl were obtained fr om the Weizmann Inter active 
Supernova data REPository (lYaron fe Gal-Yamll2012l ; WIS- 
eREPfl The best match is the +97 d spectrum of SN 1998S, 
which is quite similar to the +124 d spectrum of SN 2009ip. 
However, Ca is much stronger in SN 1998S, while the He 
lines are weaker compared to SN 2009ip. 

3.2 NIR spectroscopy 

Aside from the optical spectroscopy, we obtained a sequence 
of NIR spectra of SN 2009ip as shown in Fig. [7] These spec- 
tra are dominated by H, with clear detections of the Paschen 
lines and Br7 at all epochs. We also see strong He emission 
at 1.08u.m, and more weakly at 2.06um. Besides these el- 
ements, there are no other lines present. Most notably, we 
see no signs of O in the J-band region, while Mg at 1.5u.m 
is similarly absent. The spectra do not show any change in 
their slope over the period between October 16 and Decem- 
ber 20, in agreement with the flat J — H and H — K colours 
seen in Fig. 1151 

Our highest resolution NIR spectrum is taken with 
LBT+LUCIFER on October 25, and is shown in detail in 
Fig. Pa/3 and 7 have a FWHM of 330 kms -1 , slightly 
higher than the narrow component of the Balmer lines in 
the optical. Interestingly, we see structure in to the blue 
wing of the Pa lines, with emission components at veloci- 
ties of -1950 kms -1 and -1550 kms -1 relative to Pa7 and 5 
respectively. This structure, which is not visible in the pre- 
ceding or following SOFI spectra, is contemporaneous with 
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Figure 5. Comparison of the September 2009 spectrum of SN 2009ip presented by P13, and the final spectrum we obtained in December 
2012. The same lines are present at both epochs, albeit with higher velocities in 2012. In general, the spectra are quite similar, although 
we note that the Ca and [Fe] emission is more pronounced in 2012. 
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Figure 7. The sequence of NIR spectra of SN 2009ip. Regions of strong telluric absorption where the S/N of the spectrum is poor are 
indicated with a © symbol, the rest wavelengths of the H and He lines visible in the spectra are also indicated. 
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Figure 6. Optical spectra of SN 2009ip at represent ative epochs, 
compa r ed to a sampl e of Ty p e Iln SNe. Dat a are fr om| Fassia et al.l 
1120011 ), IZhane et at] <2012h , llnserra et all (|2012| V The phase of 
each spectrum is indicated in the plot; for the spectra of SN 2009ip 
the phase from the start of the 2012a event (Aug 1) is given, 
together with the phase from the start of the 2012b event in 
parentheses. 



the emergence of the blue shoulder in Ha in the optical. We 
note however, that the blue shoulder at in the optical spectra 
was present at velocities of ~3000 kms -1 , higher than what 
is seen in the NIR. We also find evidence for a weak P-Cygni 
absorption in Pa5 and possibly Pa/3, with minima of -3000 
kms -1 for the former, but only -900 kms -1 for the latter. 
The minimum in Pa/3 is much lower than that observed in 
other lines, and is comparable to the typical velocities seen 
in LBV outbursts. We must caution however, that the S/N 
is not as good here as at shorter wavelengths, and so the 
identification of this component is more tenuous. 

We see a broad absorption at 1.043u., which has a ve- 
locity minimum of 14500 kms -1 with respect to Pa7. We do 
not see a corresponding minimum in Pa/3 or S, and so con- 
sider it more likely that this is associated with the strong 
He I line at 1.083 u, in which case the absorption minimum 
lies at 11500 kms -1 . This is similar to what we measure in 
the Na/He blend in our optical spectra from October 30. 



9000 10000 11000 12000 13000 
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Figure 8. Close up of the LUCIFER spectrum described in 
the text, with the emission components in the blue wing of the 
Paschen lines indicated. 

3.3 The line profiles of SN 2009ip 

The H and He lines in SN 2009ip are characterized by narrow 
emission on a broad base, and multiple broad absorptions 
at velocities on order of 10 4 kms -1 . Below we discuss each 
component of the Balmer lines in turn. 

• Narrow H emission 

Similar to what is observed in an H n region, the narrow 
lines arise from the CSM. While the H in an H n region 
is ionized from massive OB stars, in the case of Type Iln 
SNe the CSM is ionized by either the flux of the SN and/or 
the emission of the cool dense shell. The velocities in the 
narrow lines (~500 kms - 1 ) are similar to th at seen in many 
Type Iln SNe and LBVs (|Kiewe et alj|2012ft . although lower 
than those see n in some extreme ex amples of the class of 
SN impostors l|Pastorello et al.ll2010l ). The narrow lines are 
present at all epochs in SN 2009ip, although they are de- 
creasing in strength with respect to the broad component 
over the period of the 2012b event. We see no narrow P- 
Cygni absorption in Ha from a wind, as shown in Fig. 1101 
T his is in contrast to the sample of interacting SNe studied 
bv lKiewe et all (|2012T ). 

We measured the central wavelength of the narrow com- 
ponent in Ha in all our spectra by fitting a Gaussian to the 
line, taking care to exclude the broad base from the fit. A 
plot of central wavelength as a function of time is shown in 
Fig. [5] which reveals a small but significant trend towards 
redder wavelengths with time. 

The slight shift of the narrow component of Ha to redder 
wavelengths is unusual, although we note that a similar ef- 
fect was seen bv lKankare et al.l (|2012T ) in SN 2009kn. A shift 
to the blue could be easily explained by the formation of dust 
in the ejecta. While one may invoke dust destruction to ex- 
plain the opposite effect, the shift to the red is quite gradual, 
and does not appear to coincide with the peak of the UV 
flux which might be expected to destroy the dust. Further- 
more if there was also significant dust being destroyed, a 
NIR infrared excess which disappears over time would also 
be expected, although we do not observe this in our photom- 
etry. Dust formation in a clump of material moving towards 
the observer could be an alternative explanation, although 
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Figure 9. Wavelength shift of the narrow components of Ha and 
H/3, relative to the rest frame of NGC 7259. Filled points indicate 
measurements obtained from the EFOSC2 spectra, which were 
pipeline reduced in a homogeneous fashion and had their wave- 
length calibration checked against night sky lines. Open symbols 
indicate measurements made on spectra from TNG+DOLORES 
and the ANU 2.3m+WiFeS. 




-1500 -1000 -500 500 1000 1500 

Velocity (km s~ ) 

Figure 10. Ha from our highest resolution spectrum, obtained 
on Oct 23. We see no evidence for a narrow (< 250 kms ) P- 
Cygni absorption at high resolution. 

we do not see a strong NIR excess in our photometry which 
may be associated with dust formation. 

• Broad emission The broad emission is approximately 
centered on the rest wavelength of Ha. While broad emission 
is usuall y associated w ith emission from ejecta at a large 
velocity, IChugail l|2001f ) first suggested for the Type Iln SN 
1998S that an alternative explanation could be scattering o f 
photons; a similar result was found bv lDessart et all ( 2009) 
for SN 1994W. Repeated Thomson scattering of photons by 
electrons in the dense CSM can give rise to broad lines which 
have a Lorentzian profile, which can be seen most clearly in 
the Oct. 8 spectrum of SN 2009ip. An alternative mechanism 
to create broad lines is boxy emission from interaction, as 
noted bv lChevalier fc Franssonl (|l994f ). 

The broad emission becomes notably asymmetric after 
November 4, approximately coinciding with the first bump 
in the lightcurve. This raises the possibility that the bump 



is due to the ejecta encountering a denser region of CSM ex- 
pelled by the progenitor star. If the density is not uniformly 
enhanced in all directions, then we may expect asymmetries 
to arise in the profile of Ha. 

• Broad absorption. We see broad absorption extend- 
ing out to a velocity of -16,000 kms -1 , with a flat bottomed 
profile. While a shallow, flat-bottomed absorption trough is 
unusual in SN and LBVs, it can be expl ained by the pres- 
ence of a thin shell o f absorbing mater ial l|jefferv fc Branch! 
1 1990( 1 Interestingly, Ijefferv fc Branch! find that in such an 
arrangement the maximum absorption in the line is 20 per 
cent of the flux in the continuum, which is in good agree- 
ment with what we see in the Ha trough in the EFOSC2 
spectrum from November 14. 

The presence of the broad Ha absorption in the spectrum 
at peak (Fig. [2| indicates that the optical depth of the cool 
dense shell cannot be much greater than a few. Alternatively, 
we require a clumpy CSM-ejecta interaction, where the high 
velocity ejecta is visible through windows in the cool dense 
shell, and which could be caused by instabilities in the ejecta 
l|Chevalier et alJll992h . 

A final interesting point of note in SN 2009ip is that the 
broad absorption appears stronger in H/3 than Ha. This is 
in contrast to what is expected, given that the optical depth 
in Ha is always higher than in H/3. A possible explanation 
for this could be broad Ha emission "filling in" some of the 
absorption in Ha. At high densities, H/3 is quenched as the 
n=4 level is depopulated to n=2 via Paa. In this case, one 
would also expect to see strong Paa emission, but as this 
line is in the telluric band at 1.87u., it is difficult to test this 
hypothesis. 



4 PHOTOMETRY 

4.1 Lightcurve evolution 

The photometry of the 2012a event, together with the initial 
rise of the 2012b event have been presented by P13. There 
is some overlap between these data, and our first photome- 
try obtained from NTT+EFOSC2. We find excellent agree- 
ment with the magnitudes reported by P13 in the region of 
overlap (Fig. 1110 ; both sets of photometry agree to within 
their uncertainties, and show no systematic offsets. We have 
als o compared our l ightcu rve to the BRI data presented 
by iMauerhan et al.l (|2013T ), and find in general fair agree- 
ment between the two data sets within the uncertainties, 
and again see no systematic offset. The region of the bump 
in the lightcurve (arou nd Nov 7) is the only discrepancy, 
where IMauerhan et al.l find a similar increase in 7, a larger 
increase in B, and no increase in R, although the reason for 
this disagreement is unclear. 

The lightcurve in all filters exhibits a decline from the 
bolometric maximum of the 2012b event (around October 
3, MJD 56203), until the end of the observing season on De- 
cember 29 (MJD 56290) when SN 2009ip disappeared be- 
hind the Sun. The decline is the strongest in the UV bands, 
with a drop of 8 mag from peak in uvw2 over 50 days, com- 
pared to only a 2 mag drop over the same period in /. During 
our last observation of SN 2009ip, taken on December 29 in 
V, the transient had faded to My=-13.1, 5 magnitudes be- 
low peak. Even at this phase the SN is still 4.5 mag brighter 



© 0000 RAS, MNRAS 000, 000-000 



SN 2009ip 13 



Date 

Aug 01 Aug 21 Sep 10 Sep 30 Oct 20 Nov 09 Nov 29 Dec 19 




56140 56160 56180 56200 56220 56240 56260 



MJD 

Figure 11. Optical (square points), UV (round points) and NIR (triangle points) photometry of SN 2009ip for the 2012a and 2012b 
events. Open points are taken from the data pub l ished by P13, while filled points are new data presented in this paper. We also include 
the JHK mag nitudes reported for Oct 29 by[lm| ||2012T ). Also indicated are the quiescent F606W (R) magnitude of the progenitor from 
I Foley et al. ( 2011) (solid arrow, offset by -2 mag to match the R-band magnitudes of SN 2009ip), and the expected 0.98 mag per 100 
days decline rate in V given by the radioactive decay of 56 Co (dashed line). 

in R t han the quiescent magnitude reported bv lFolev et al.l 
l|201ll) . F606W=21.84±0.17. 

Around November 1, we see a clear bump in the 
lightcurve of SN 2009ip, where the transient brightens by 
0.2-0.5 mag in all filters. Curiously, the bump appears to de- 
velop first in the redder filters. We show a zoomed in section 
of the lightcurve of SN 2009ip in Fig. 1121 where we have fit- 
ted natural splines to the photometric points between MJD 
56230 and 56245 in each filter (taking the points immedi- 
ately prior and subsequent to this range as the boundary 
values, and weighting all points equally). From the spline 
fits it appears that the bump becomes visible first in R and 
7, and lags by a few days in the Swift UV filters. The bump 
also appears to be more pronounced in the optical; while we 
have probably missed the peak in RI, UBV show a clear 
rise, as opposed to uvwl,uvm2,uvwl which merely flatten 
for a few days. 

A second bump in the lightcurve can be seen at the 
start of December (MJD ~56265). This bump is less pro- 
nounced than the one seen around November 1, with the op- 
tical bands merely flattening for a period of around a week, 
rather than increasing in brightness. Unfortunately the sec- 
ond bump coincides with the point at which SN 2009ip is 




56230 56235 56240 56245 

MJD 

Figure 12. The region of the first bump in the lightcurve of SN 
2009ip (around Nov 1), along with spline fits to the data. 

too faint to observe in the UV with Swift, while the R and 
/ bands are sparsely sampled to clearly identify the peak of 
the bump. The December bump lasts for a similar duration 
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Figure 13. Comparison of the absolute R-band lightcurve of 
SN 2009ip with a selection of interacting transients (SN 2008S, 
t =54485.5, A fl =1.6, /i=28.7, iBotticelfa et al I|200S| ; SN 1998S, 
t =50874.7, A fl =0.56, u=31.15. iFassia et al.l l2000l: SN 2010jl, 
t maa; =54887.5, A H =0.06, u=33.44 jZhang et alj2012t SN 200 5ip, 



t disc =53679.2, A fl =0.011 
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2009kn, t =55115.5, A R -- 
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34.23, lKankare et alj|2012l) 
Lightcurves are plotted with respect to the explosion epoch (to), 
or where a good estimate of this is unavailable, to the discovery 
date (tdisc) or epoch of maximum light (t max ). Distances and 
extinctions are taken from the literature. 



to the one in November - perhaps a week or less - and is fol- 
lowed by a continuing decline in the lightcurve afterwards. 
From around December 17 (MJD 56278) the lightcurve ap- 
pears to flatten until the last epoch on December 29. 

In Fig. [13] we compare the .R-band lightcurve of SN 
2009ip over August - December 2012 to a sample of inter- 
acting transients taken from the literature. The peak mag- 
nitude of SN 2009ip, ~ —18, is comparable to that of SNe 
2005ip and 2009kn, but much fainter than that of Type Iln 
SNe such as SN 2010J1 and 1998S. SN 2009ip also has the 
most rapid decline from peak of any objects in the sample. 

The colour evolution of SN 2009ip is presented in Fig. 
[Til For the 2012a event, and up to the peak of the 2012b 
event, the colour in all optical bands is flat. As the SN fades 
from the peak of the 2012b event the colour becomes redder 
in all bands as the photosphere cools. The colour evolution 
from photometry is in accord which that expected from the 
spectra, which show a similar slope until peak. The B — V 
colour shows a bump around the same time as the Nov 7 
undulation in the lightcurve, which is to be expected given 
the time difference between the bump in different filters, as 
discussed in Section 14.21 A comparison to other interact- 
ing transients (SNe 1998S and 2009kn) shows a good match 
with SN 2009ip, with the only notable difference being a 
bluer U-B color for the latter compared to SN 2009kn. The 
NIR colours of SN 2009ip are shown in Fig. [15] As in the 
optical, we see a similar color evolution to SNe 2009kn and 
1998S. The evolution of SN 2009ip is consistent with a cool- 
ing blackbody, and shows no evidence for a NIR excess which 
would signal dust formation. 

We constructed a bolometric lightcurve for SN 2009ip 
from our UV, optical and NIR photometry. The measured 
magnitudes were converted into fluxes at the effective wave- 
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Figure 14. Optical colour curves for SN 2009ip during the 2012a 
and 2012b events. Open symbols are from P13, w hile closed sym- 
bols are from this pape r. SNe 2009kn and 1998S jKankare et al] 
120121; I Fassia et al]|20o"oT ) are also shown, these objects have been 
shifted so that their explosion epoch is at MJD=56140. 
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Figure 15. NIR colour curves for SN 2009ip during the 2012a 
and 2012b events. Open symbols are from P13, while closed sym- 
bols arc from this pap er. We al so include the JHK magnitudes 
reported for Oct 29 bv llinl <2012h. SNe 2009kn and 1998S are als o 
shown for comparison llKankare et al] 120121 ; I Fassia et al] 12000?) , 
shifted so that their explosion epoch would be at MJD=56140. 
All colour curves have been corrected for reddening. 



lengths of the filter used and corrected for extinction, be- 
fore the resulting spectral energy distribution was integrated 
over the wavelengths covered with a linear interpolation, 
assuming zero flux beyond the limits of the integration. 
Fluxes were computed at epochs where there was an ob- 
servation in V; where observations in other bands were not 
taken contemporaneously, their magnitudes were estimated 
by interpolating the light curves using low-order polynomi- 
als, or were extrapolated assuming constant colours. The 
bolometric flux at each epoch was converted to a luminos- 
ity given the adopted distance to SN 2009ip. The resulting 
bolometric lightcurve is shown in Fig. 1161 The bolometric 
lightcurve of the 2012a event is dominated by optical light. 
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Figure 16. The uvoir bolometric lightcurve of SN 2009ip (black), 
together with the contributions of the optical, UV and NIR pho- 
tometry 
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Figure 17. The uvoir bolometric lightcurve of SN 2009ip (black), 
together with two oth er well observed Type Iln SNe ; SN 2 009kn 
llKankare et al.l I2012T) and SN 2010jl dZhang et alj|2012h . The 
lightcurves of SNe 2009kn and 2010jl are relative to their ex- 
plosion date, while that of SN 2009ip is relative to the start of 
the 2012a event. 



For the rising phase of the 2012b event however, the con- 
tribution to the lightcurve from the UV is larger than from 
optical wavelengths. The NIR contributes very little flux to 
the lightcurve of SN 2009ip for almost the entire duration 
of the 2012a and 2012b outbursts, amounting to only 20 
per cent of the total at the time of our last epoch of NIR 
photometry. 



4.2 The undulations in the lightcurve of SN 
2009ip 

The presence of bumps in the lightcurve of SN 2009ip is un- 
usual (but not unpreced ented, as they were al so noted in 
the Type Iln SN 20051a; IPastorello et alj|2008h . To repro- 
duce these short bumps with the interaction of a few solar 
masses of ejecta (as might be ejected in a Type Iln SN) with 
the CSM may be difficult, as a massive SN ejecta should be 



more homogeneous and not give rise to small scale variations 
in the lightcurve. It is also possible that the bumps are peri- 
odic, although monitoring over a longer timescale is needed 
to confirm this. In this case, the bumps could possibly be 
caused by spiral density structures resulting from a passage 
of a binary companion. 

Another surprising aspect of the bumps is that they 
appear to peak first in the redder bands. Dust is unlikely to 
be the cause of the lag in the UV, as while it may cause the 
bump to be less pronounced at bluer wavelengths, there is no 
obvious reason why the peaks of the lightcurve in different 
filters would not coincide. 

An alternative explanation for the shift in peak with 
wavelength is successive scattering by metal lines in the UV. 
The peak in U and the Swift UV filters is approximate 1- 
2 days after the peak in BVRI. To explain this delay with 
an additional travel time for UV photons implies that they 
cover a distance of 170-340 AU more than the optical and 
NIR photons. For a characteristic radius of ~40 AU to the 
cool dense shell, this requires each photon undergoes multi- 
ple scattering events. It is also important that the photons 
(or at least a high percentage of them) scatter from the lines 
(which will give rise to re-emission at the same wavelength), 
rather than causing the atom to fluoresce (which would shift 
flux towards redder wavelengths) or undergo collisional de- 
excitation. 

There is no marked change in the spectra of SN 2009ip 
around the time of the first bump. The [Ca n] beings to 
emerge at this epoch, which may indicate a lower density, 
but it is difficult to envisage a decrease in density causing a 
rise in luminosity. 

4.3 Is there Ni in SN 2009ip? 

The final epochs of photometry obtained for SN 2009ip be- 
fore it disappeared behind the Sun appear to show a flatten- 
ing in the redder bands. Given that the lightcurve displayed 
undulations prior to this (as discussed in Section f4.2p . it is 
entirely possible that the apparent flattening is merely the 
start of a third bump in the lightcurve, rather than the be- 
ginning of any 56 Ni-powered tail phase. We note that the 
overall decline in the lightcurve from peak (4.5 mag over 80 
days in V) is much faster than that expected from the decay 
of 56 Co (0.98 mag per hundred days). 

From our final photometric point on December 29, we 
can estimate an upper limit to the mass of 56 Ni synthesized 
in SN 2009ip (assuming it was indeed a genuine SN). We 
measure a bolometric luminosity of log(L)= 40.8 erg s -1 fo r 
SN 2009ip on December 29. Using eq. 2 of lHamu 



rg is tor 

v (2003), 



and assuming that SN 2009ip exploded on 2012 Aug 1, we 
calculate a 56 Ni mass of <0.02 Mg. This is an upper limit to 
the mass of 56 Ni in SN 2009ip, as it is based on the earlies t 
plausible date for a SN explosion (jMauerhan et all 120131 ') . 
If the 2012a was not a SN, but rather the star exploded 
at the start of the 2012b event, t hen the 5 6 Ni m ass must 
be less than O.O1M0. Furthermore, iHamuvl (|2003T ) assumes 
that the sole energy source powering the luminosity of the 
tail is radioactive decay, whereas in the case of SN 2009ip 
there is at least some luminosity arising from interaction (as 
evidenced by the persistence of narrow emission lines until 
this epoch). 

One caveat to this calculation is that the forma- 
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tion of dust can serve to mask the late time luminos- 
ity of SNe at optical wavelengths, and cause their 56 Ni 
mass to be underestimated. However, in the case of SN 
2009ip, we see no evidence for dust formation in either 
spectroscopy (no apparent blue shift of lines) or in pho- 
tometry (no strong NIR excess). A second caveat is that 
we assume that 7-rays from the decay of 56 Co are fully 
trapped within the e j ecta. For the well-observed SN 1987A 
jSuntzeff fc Bouchetl I1990T ). approximately 2/3 of 7-rays 
were still trapped at 500 d post-explosion. The ratio of the 
trapping timescale in SN 2009ip to that in SN 1987A will 
scale as (Mo9i p / 'M$r a) ' 5 (uo9ip/ u 87A) _1 . If we assume that 
the velocities are roughly similar in the two events, and that 
the ejecta mass of SN 1987A was 12Mq, then this scaling 
implies that for 7-ray leakage to be significant at 100 or 150 
d (depending on assumed explosion epoch), the ejecta mass 
in SN 2009ip should be less than 0.5 or 1.1 Mq respectively. 
Such a small ejecta mass is much lower than that typically 
seen in SNe, and appears inconsistent with the core-collapse 
of a massive star. Even with these caveats, the restrictive 
upper limit on the mass of any 56 Ni in SN 2009ip appears 
robust. 

4.4 The energetics of SN 2009ip 

Integrating under the uvoir bolometric lightcurves of the 
2012a and 2012b events, we find an observed radiated energy 
of 7.8xl0 47 ergs for the former and 1.8xl0 49 ergs for the 
latter. In a normal CCSN, the radiated energy over the first 
one hundred days is typically comparable to that radiated 
in the 2012b event, ~ 10 49 ergs. 

While the large energy of the 2012b event is suggestive 
of a genuine SN explosion, with core-collapse as the energy 
source, it is interesting to consider the possibility that the 
high luminosity is in fact completely due to the efficient 
conversion of kinetic energy from colliding ejecta. From the 
broad absorptions in the Balmer lines we see material at 
velocities up to 12,000 kms -1 , however we adopt 7,500 as 
a more typical velocity for the bulk of the ejecta. If we as- 
sume 100 per cent efficiency in converting kinetic energy to 
radiation, then we only require 0.05 M0of material at 7500 
kms -1 to give rise to the luminosity of the 2012b event. 

To check whether the assumption of high efficiencies 
in the conversion of kinetic energy is correct, we apply the 
Rankine-Hugoniot conditions to the shock. For ejecta mov- 
ing at v e j hitting a reverse shock at v ra , the efficiency e in 
converting kinetic energy to luminosity is given by 

vl] ~ (VVej + (1 -V)Vrs) 2 

6 ^ V 2 - { ' 

ej 

where £ is the porosity {Aread, umj>3 / AreacsM) to ac- 
count for dumpiness in the CSM, and rj is the ratio of post- 
and pre-shock velocities in the rest frame of the shock. 77 de- 
pends on the equation of state of the gas, ranging between 
1/4 and 1/7 for a perfect gas and a radiation dominated gas 
respectively. If we take the reverse shock to be moving at 
approximately the velocity of the photosphere derived from 
our blackbody modelling, 1200 kms -1 , and the ejecta to be 
moving at a conservative estimate of 7500 kms -1 , then we 
derive efficiencies of around 90 per cent. These efficiencies 
will be reduced by £, although even if 50 percent of the CSM 



shell was clumps and 50 percent were holes, then we would 
still have e ~ 0.45. In this case, we would still only require 
~0.1 MQof ejecta to power the interaction luminosity. 



5 MODELLING THE CSM AROUND SN 2009IP 

As a first attempt at understanding the physical condi- 
tions in SN 2009ip, we have constructed a spherically sym- 
metric CSM model of SN 2009ip with a density profile 
R -2 as detailed in Appendix fA] While the real structure 
of SN 2009ip is undoubtedly much more complex, with 
multi dimensional effects and prob ably multiple shocked 
regions Dwark adas fc Gruszkd (|2012l ). less complex models 
can hopefully be used to derive order of magnitude results. 

We perform a least-squares fit for black-body radii and 
temperatures (4TvR 2 TrB\(T t ) = L\) to the observed, de- 
reddened, UBVRI photometry for selected epochs (typi- 
cally the ones with full photometric coverage). The fit for 
the epoch we model in detail in Section [S] (Oct 8 2012) is 
shown in Fig. [TH] as an example, while the time evolution 
of the derived radius and temperature is shown in Fig. 1191 
The radius at the start of the 2012b event derived from 
a black body fit is also consistent with the maximum ra- 
dius derived from light travel time con siderations using th e 
high cadence photometry presented by iPrieto et al.l l|2013h . 
The evolution observed over the 2012b event is quite sim- 
ilar to that seen in SN 1994W, which peaked at 15000 K 
and declined to 7000 K over 70 d, but hott er than that of 
SN 2 009kn, which only peaked at 10000 K (|Kankare et all 
l2012j l However, as the first spectra of SN 2009kn were ob- 
tained later than for SN 1994W, it is possible that the 
maximum temperature was in fact higher. Comparing to 
the radius evolution of SNe 1994W and 2009kn, we see 
that the initial photospheric radii of both these objects 
(R~8,500R Q ) was larger than that of SN 2009ip during the 
2012a event (R~3000Rq), but that they expanded with sim- 
ilar velocities to the 1200 kms -1 observed in the 2012b 
event (vsni994w =1000 kms -1 , usni994w =1000 kms -1 ; 
«SN2009kn =600 kms -1 ). 

We use the photospheric temperatures and radii derived 
as radiative boundary conditions for computing the ioniza- 
tion and excitation balance of hydrogen in the CSM as de- 
scribed in Appendix fA] We do not solve for temperature, 
but assume T = T*, which should be a reasonable approxi- 
mation. 

In particular, we investigate the first epoch after maxi- 
mum light when we have an optical spectrum, MJD = 56209 
days, when T, « 12,600 K and R* « 44 AU. The lack of 
any strong broad lines from the ejecta at this epoch, com- 
bined with the large photospheric radius is possibly due to 
the formation of a cool dense shell at the interface between 
the outer ejecta and CSM, which becomes optically thick in 
the Paschen co ntinuum for a few weeks if the CSM density 
is high enough (|Chugaill200l1) . Outside this will be a region 
of shocked CSM, and outside the forward shock a region 
of unshocked CSM, radiating due to ionization by the pho- 
tospheric radiation field (as well as from X-rays from the 
forward shock, which we however do not include here). An- 
other possible explanation for the lack of broad absorptions 
is that the fast material is efficiently ionized by the emission 
from the shock region. 
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Figure 18. The (dereddened) spectrum at Oct 8 2012 (black, 
solid), together with (dereddened) photometry (red points) and 
black-body fit to the photometry (blue, dashed). The derived tem- 
perature is 12650 K. 
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Figure 20. Ratio of Ha and photospheric luminosity, on Oct 8 
2012 (first spectrum after maximum light). The observed value 
(3 X 10 — 3 ) is marked with a dashed line. The models with a 
velocity-gradient are terminated at the density where the Sobolev 
limit is no longer reached (at high densities they behave as the 
model with zero gradient). 
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Figure 19. Upper panel: The temperature evolution of SN 2009ip 
as measured from black-body fits to optical photometry at se- 
lected epochs. The approximate dates of the 2012a and 2012b 
events are indicated. Lower panel: The photospheric radius of SN 
2009ip as determined from blackbody fits to optical photometry. 
The solid line is a least squares fit to the radii measured for the 
2012b event only, while the dashed line is the extrapolation of 
that fit to the earlier 2012a event. The slope of the fit is 1200 
kms . 



We explored the sensitivity of the model to having a 
region between the photosphere and the CSM, but do not 
find a significant difference. 

For the inner parts of the CSM, the photospheric ra- 
diation field is strong enough that radiative processes are 
important not only for the ionization balance but also for 
the NLTE solutions to the excitation balance. 

We may therefore not use standard recombination the- 
ory, which assumes that the external radiation field is ne- 
gieglible for the statistical equilibrium (as it is in normal HII 
regions where the dilution factor (R/R*) 2 is large). 

We explore models with different assumptions for the 
velocity gradient in the CSM (as discussed in Appendix [A"| . 
Fig. [20] shows that a model with a zero velocity gradient 



(i.e. a constant velocity flow throughout the emitting region) 
is not successful in reproducing the observed strength of 
This is a consequence of the efficient line trapping 
in a constant velocity flow; when the density becomes high 
enough to trap any significant fraction of the photospheric 
radiation (n > 10 9 cm -3 ), Ha has such a high optical depth 
that n — 3 (the level creating Ha) becomes depopulated 
by collisions and photoionization rather than spontaneous 
decay. 

Models with velocity gradients dV/dr > dV/dro (~ 
1.4 x 10 -7 ; our estimated velocity gradient, found by di- 
viding the CSM velocity by the scale size of the region), 
produce more promising values of the Ha luminosity for den- 



sities n ~ 10 y - 10 iu cm~ 3 . Turning to Fig. 21(a) and 21(b) 



we see that these models also produce correct line ratios of 
L(Ha)/L(Hb) « 2 and L(H<y)/ L(H /3) w 0.45. As Fig.lHTc) 



shows that Ha is marginally optically thick in these models, 
the reason for the damped Balmer decrement seems to be 
that Ha has an optical dept h of a few, but H/3 is optically 
thin, a scenario discussed bv lLevesque et all (|2012l ). 

Fig. 21(d) shows that all models predict electron scat- 
tering optical depths of order unity for densities between 
n = 10 9 - 10 11 cm . Since the Balmer lines show distinct 
broad wings that likely arise du e to electron scattering with 
optical depth of order a few llChugail l200ll; IChugai et al.l 
12004 iDessart et ai1l2009l ; ISmith et al.ll2010bl ). this further 
strengthens the idea that the CSM density lies in this inter- 
val. 

As made clear by the previous figures, the high density 
n > 10 13 cm -3 proposed by Levesque as one possibility for 
the low Balmer decrement is not succesful in reproducing the 



7 The observed value was obtained by measuring the relative 
amount of flux in the line compared to the total flux in the spec- 
trum (3300 - 10000 A), plus estimated fluxes outside the spectral 
range from the black-body fit. This gives a fraction of 5 X 10 — 3 . 
Assuming no flux outside the spectrum, the value is 0.012. 
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(a) Hct/H/3 ratios as a function of density in our models for the 
CSM lines around the maximum of the 2012b event. The observed 
value is indicated with a dashed line. 
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(b) Same as Fig. |21(a)| but for the H7/H/3 ratio. 
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(c) Model Ha (solid) and H/3 (dashed) optical depths around the (d) Electron scattering optical depths in our models for the CSM 
peak of the 2012b event. lines around the peak of the 2012b event. 



Figure 21. Results from our modeling as described in Section[5] 



observables. At these high densities the Balmer lines have 
high optical depth for any reasonable assumption about the 
velocity field, and as the level populations also approach 
LTE, the line ratios approach their blackbody values. As an 
example, the H^y/H/3 ratio goes towards 1.3 at high densities 
(and T — 12000 K), whereas the observed value is ~ 0.45, 



as reproduced by the NLTE models around n = 10 



10 1 



cm . A high-density model would need T ~ 2500 K to re- 
produce Hy/H/3 ratio, but at that temperature the Balmer 
decrement has value 4.8, much too high. The other model re- 
sult that clearly rules out a n > 10 13 cm -3 scenario is that 
that the electron scattering optical depth would be much 
smaller than unity. 

The number density for a constant mass loss rate at a 
radius r is 



»(r) = 



M 



Am p 4nr 2 v(r) 



r ln io -3 
= 5-10 cm 



(40 Au) 



M 



1 M Q yr-i 



V 



10 3 km s- 1 



(2) 



where M is the mass-loss rate, A is the mean atomic 
weight (1.6 for a normal H/He mix), m p is the proton 
mass, and v(r) is the wind velocity. The derived density 
of n ~ 10 9 — 10 10 therefore suggests as mass-loss rate of 
M ~ 0.01 — O.IMq /year. We caution however, that this is 
likely only correct to an order of magnitude, due to the un- 
certain assumption of a constant mass-loss rate. 

Using more sophisticated modelling, iDessart et alj 

(|2009h find a similar Balmer decrement of ~2 for SN 1994W 
from a model with a de nsity of 10 9 — 10 10 cm " 3 and with a 
velocity of ~ 800km s ~ 1 , [ChugaT et all (|2004T ) also require a 
similar density in order to reproduce the electron scattering 
wings. That these values are in close agreement with our 
simpler calculations lends credence to our approach. 



6 ANALYSIS 

There are thee questions which are key to understanding 
SN 2009ip. Firstly, what is the significance of the surpris- 
ingly remote location of SN 200ip; secondly, were either of 
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Figure 23. Enclosed R flux from NGC 7259 as a function of ra- 
dius, normalized to total flux. The arrow indicates the position of 
SN 2009ip, which lies outside the detectable extent of the galaxy. 



Figure 22. BVr color composite of SN 2009ip (marked with cross 
hairs) and NGC 7259 obtained with the Liverpool Telescope + 
RATCAM on 20f2 November 13. Scale and orientation are as 
indicated. The 0.5' scale bar corresponds to a physical distance 
of approximately 3 kpc at the distance of NGC 7259. 



the 2012a/b events a SN explosion or a mega-outburst of 
an impostor; and thirdly, what powers the luminosity of the 
2012a/b events? As these questions are linked, we will ad- 
dress them in turn. 



6.1 The lonely fate of SN 2009ip 

One of the most striking features of SN 2009ip is its iso- 
lated location, as illustrated in Fig. 1221 SN 2009ip is at a 
projected distance of ~0.75' from the center of NGC 7259, 
placing it at a deprojected physical distance of of 4.3 kpc. To 
measure the radial extent of the galaxy, we took the deepest 
R band images of SN 2009ip taken with EFOSC2 under 1" 
seeing or better (from 2011 June 23), and stacked them to 
create a deep 1200s exposure. No faint spiral arm or other 
obvious structure can be seen at the SN location, in fact the 
nearest extended source to SN 2009ip in NGC 7259 is ~10 
" distant. We measure a surface brightness at the SN loca- 
tion of R=25 mag arcsec 2 . The radial extent of NGC 7259 
has been measured at 0.6' to 0.78' for surface brightness of 
25 and 27 mag arcsec 2 resp ectively (as listed in NED, from 
lLauberts fc Valentiin|[l989l ). Hence although there is little 
star formation in the vicinity of SN 2009ip, it appears to 
reside within the optical disk of the galaxy. 

After masking the nearby galaxy 2MASXi J2223042- 



285647 in our deep image, we performed aperture photom- 
etry on NGC 7259 using circular apertures between 5 and 
250 pixels (i.e. to a maximum radius of 1'). The sky back- 
ground was measured in a blank region of the image sev- 
eral arc minutes distant from NGC 7259 to avoid any un- 
intended contamination. We plot the enclosed flux from 
NGC 7259 as a function of rad i us in Fig. [23] in a simi- 
lar fashion to lAnderson fe Jame 1 j2009l ). The fraction of R 
band light within the position of SN 2009ip, F rR , is ~1. 
lAnderson fc James! found that the mean F rR for Type Iln 
SNe and SN impostors were slightly larger than for Type 
IIP SNe, and that there was a deficit of both impostors and 
Type Iln SNe in the central regions of galaxies. While SN 
2009ip is in agreement with this trend, it is an extreme ex- 
ample, and indeed has a larger value of F rR than any of 
the interact ing transients in their sample. This is also in 
accord with I Anderson et al.l ((2012) , who found that SN im- 
postors and Type Iln SNe are the least correlated with on- 
going star formation as measured from Ha and UV imaging. 
It is w orth noting, how e ver, t hat the "SN impostor" cate- 
gory of lAnderson et all (|2012l ) is extremely heterogeneous, 
and co ntains quite dissimila r transients (for example SN 
2008 S iBotticella et al.ll2009l and SN 2002kg, iMaund et all 
2006:; IWeis fc Bomans!l2005l ). Nonetheless, the location of SN 
2009ip is quite unique among bo th Type Iln SNe, and SN 
impostors. [Anderson et al.l l|2012l ) have used the argument 
that the lack of correlation of Type Iln SNe with Ha emis- 
sion implies they have lower mass progenitors than the bulk 
of the CCSN population. While SN 2009ip is certainly iso- 
lated from any strongly starforming regions, it has (without 
doubt) a high mass progenitor star. 

The HST+W FPC2 F606W images presented by 
iFolev et all (|201lh are also of use for con straining the en - 
vironment of SN 2009ip. Using hstphot (|Dolphin I [20001 ). 
we detect only a single sourcqj within a 4" radius of SN 
2009ip. We calculate a limiting magnitude of F606W>25 
mag, implying that there is only a single source with a R 



8 We are excluding here the foreground star (R~18.4) to the NE 
of SN 2009ip. 
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magnitude brighter than -6.6 within a 500 pc radius of SN 
2009ip. 

The apparently isolated position of SN 2009ip is of 
interest for two reasons. Firstly, as a LBV with a ZAMS 
mass of ~50-60 Mq, the progenitor of SN 2009ip could be 
an example of isolated starformation. While there are ex- 
amples of isolat ed massive stars in the L MC (for exam- 
ple, VFTS 682, iBestenlehner et all 1201 ll ). these are still 
within a few tens of parsec of a massive star forming re- 
gion, and in the case of VFTS 682 have been suggested to 
have been dyna mically ejected from a you ng cluster with a 
top-heavy IMF (|Baneriee fc Kroupal [20121) . For SN 2009ip, 
there is no obvious nearby cluster, even with a runaway ve- 
locity of 50km s _1 , permitting it to cover a distance of 500 
pc in 10 7 years. Secondly, the extreme radial location also 
points towards a sub-solar metallicity . App lying the metal- 
licity gradient of iBoissier fc Prantzosl (|2009l ) to the position 
of SN 2009ip, and estimating the characteristic metallicity 
of NGC 7259 from its B band magnitude, we find a local 
metallicity of 12+log(0/H) = 8.5 dex. This value is signif- 
icantly lower than t he locus of SN metalliciites derived by 
IBoissier fc Prantzosl using the same technique on a sample 
of SN host galaxies. Unfortunately there is no nearby H II 
region where we could obtain a spectroscopic estimate of the 
metallicity, while the auroral [O n] and [O ill] lines are not 
detected in our spectra precluding a direct measurement of 
the electron temperature and metallicity. 

6.2 Has SN 2009ip exploded as a SN? 

Ultimately, the most definitive proof that a star has ex- 
ploded as a CCSN is the disappearance of its progenitor in 
late time images taken after the explosion has faded. While 
it should be possible to eventually test this for SN 2009ip, 
at the time of the final photometry obtained of SN 2009ip 
at the end of December, the transient was still ^4 . 5 ma g 
brighter than the progenitor detected by I Foley et alj (|201ll ). 
If the 2012b event continues its decline rate from maximum 
of ~3.5 mag in 70 d, then by the time SN 2009ip reappears 
in mid-April 2013, we may expect it to be fainter than the 
magnitude of the progenitor if it was a genuine SN. 

Another critical piece of evidence for a SN explosion 
is the detection of nucleosynthesized material in the ejecta. 
This may be either explosively produced isotopes (e.g. 56 Ni) 
or elements synthesized during the latter stages of stellar 
evolution (e.g. 16 O). In SN 2009ip we do not see any evidence 
for freshly synthesized material; all the lines present in the 
Dec 20 spectra were also present in September 2009. We 
also see none of the usual lines which are typically seen in 
the nebular spectra of CCSNe such as [O i] 6300,6364 A, 
Mg i] 4571 A and Mg I 1.5um. Against this however, we 
must note that many Type Iln SNe do not show these lines 
until quite late phases, due to the low gas densities of ~10 9 
cm' 3 required. Given the precipitous decline of SN 2009ip 
from the 2012b peak, one may expect a reasonably prompt 
transition to the nebular phase. 

Another signature of a CCSN is the presence of a ra- 
dioactively powered tail phase, which at late times follows 
the decay rate of 56 Co. From our last photometric point, 
we set a limit of <0.02 M Q of Ni in SN 2009ip. Such a low 
ejected Ni mass is hard to reconcile with the core-collapse 
of a very massive star, although a similarly low value of 



<0.023 M was set for SN 2009kn (|Kankare et all 120121 ') . 
albeit with the cav e at tha t the tail may not be Ni-powered. 
lUmeda fc Nomotol (|2008l ) find that models with a ZAMS 
mass of 60-80 M Q produce between 0.8 and 7.9 M of Ni, 
depending on explosion energy. The low upper limit to the 
ejected Ni mass could be explained if much of the Ni did 
not escape, but rather underwent fallback onto the com- 
act remnant. Such a scenario was proposed for SN 1997D 
Turatto et al.| [T998). although in the case of fallback the 
ejecta velocities are expected to be closer to 10 3 than 10 4 
kms -1 . 

We note that a low ejected Ni mass was al so derived for 
the Type Iln S N 199 4W (|Sollerman et al.ll 19981 ). In this case, 
IChugai et af] (|2004 ) claimed that the transient arose from 
the interaction of a SN with the hydrogen e nvelope of the 
star, which had been ejected 1.5 years prior. IChugai et al.l 
proposed that the ejection of the envelope could be due to 
a Ne burning flash in the O-Ne-Mg core of an 11 Mq star; 
while the low progenitor mass would also account for the 
small ejected Ni mass. The explosion of a Super Asymp- 
totic Gian t Branch (SAGB) star was also claimed for SN 
2011ht by"~ 
sibility by 



Mauerhan et al] (120121 ') . and suggested as a pos- 
Kankare et al. I~ i|2012i ) for SN 2009kn. In the case 



of SN 2009ip however, the massive progenitor is clearly in- 
consistent with such a scen ario. An alternative model for 
SN 1994W was proposed bv lDessart et al] (|2009l ). who sug- 
gested that the collision of two ejected shells with an effi- 
cient conversion of kinetic energy to luminosity could also 
power the transient. In the colliding shells model, the pro- 
genitor star is not necessarily destroyed, although the lack 
of pre-explosion imaging coupled with its crowded environ- 
ment makes this difficult to test for SN 1994W with deep, 
late-time imaging. 

Less compelling indicators of a CCSN include high ve- 
locities and luminosities. In the case of the former, the ob- 
servation by P13 of material at velocities of ~ 10 4 kms" 1 in 
2011, long before the start of the 2012a and 2012b events, 
suggests the material can be accelerated to high velocities 
without a core-collapse explosion. The bright absolute mag- 
nitude of SN 2009ip (V=-17.7) is certainly suggestive of core 
collapse, however, as we have shown in Section 14.41 the ki- 
netic energy of a relatively small amount of material at the 
velocities observed can power the lightcurve. Furthermore, 
some of the most luminous SNe known have been suggested 
to arise from outbursts whic h do not (at least im mediately) 
destroy the progenitor star (|Wooslev et alj|2007l ). 

From the available evidence - and remaining mindful 
of the truism that "absence of evidence is not evidence of 
absence" - we suggest that thus far there is no evidence 
that SN 2009ip did indeed undergo core-collapse to produce 
either the 2012a or 2012b events. 



6.3 What powers the luminosity of the 2012a and 
2012b events? 

Any scenario for SN 2009ip needs to explain the following 
key characteristics: 

• An intrinsically luminous, a most likely very massive 
progenitor star, probably with a low metallicity. 

• Multiple eruptions over at least a three year period, 
which eject material at velocities of up to 10 4 kms -1 . 
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Figure 24. Comparison of the bolometric l ightcu rve of SN 2009ip 
to the model prediction of ISoker~fe Kashi (2013). 



• A brief, 40 day eruption. 

• A subsequent more luminous eruption, with a fast rise 
time, which decays rapidly from peak. 

• An ejected 56 Ni mass <0.02 M . 

• No observed synthesized material, either from explosive 
nucleosynthesis or late-stage stellar evolution. 

There have been several suggestions as to the energy 
and n ature of the 2012a and 2012b events. iMauerhan et al.l 
(|2013h argued that the 2012a event was the explosion of SN 
2009ip as a weak CCSN with a low ejected Ni mass, while 
the 2012b event was caused by the ejecta hitting the CSM. 
P13 further suggested that if SN 2009ip was indeed a gen- 
uine SN, then the 2012b event may be the explosion of SN 
2009ip, while the 2012a event was an outburst. However, 
P13 also advanced an alternative scenario that SN 2009ip 
has not exploded as a SN, but rather that the lightcurve and 
spectra could be explained by the collision o f two massive 
shells of material ejected by the progenitor (|Dessart et al.l 
2009). Such an erupti on could be powere d by the pulsa- 
tional pair instability (|Wooslev et al.ll2007l). A final model 
for SN 2009ip was proposed bv lSoker fc Kashil (|2013h . who 
favoured a merger of a massive L BV with a bina ry compan- 
ion, in a "mergerburst" model. ISoker fc Kashil draw com- 
parison to the enigmatic transient V838 Mon, for which a 
simi lar scenario was prop osed. 

ISoker fc Kashi (|2013l ) predicted a lightcurve evolution 
for SN 2009ip from the peak of the 2012b event (the last 
epoch where published data was available to them) as shown 
in Fig. [H The overall match to the lightcurve of SN 2009ip 
is striking; the decline rate is similar, as is the amplitude of 
the undulations. However, ISoker fc Kashil also predict that 
only a small fraction of the ejecta should have velocities in 
excess of 5000 kms -1 , while in SN 2009ip we see the deepest 
absorption troughs at higher velocities than that, over three 
months after the 2012a peak. 

As discussed in Section f6. 2 1 while a weak CCSN explo- 
sion in the 2012a event cannot be entirely ruled out, we see 
no conclusive evidence in favor of it. Furthermore, while a 
weak explosion could explain the low 56 Ni mass, we would 
not expect to observe the high ejecta velocities seen. If the 
2012b event was the SN, then the upper limit to the S6 Ni 



mass is even more restrictive. Whether a ~60 M© star can 
explode, suffer sufficient fallback to hide almost all its Ni, 
but still retain a strong enough explosion to drive the ejecta 
outwards at upwards of 10 4 kms -1 remains to be seen. 

In the canonical model for the Type Iln SN, the lu- 
minosity is powered by interaction between ejecta and the 
CSM. At the start of the 2012b event, the radius of the pho- 
tosphere is around 30 AU, and is expanding with a velocity 
of 1200 kms -1 , as can be seen in Fig. 1191 If the observed pho- 
tosphere is formed at the cool dense shell, then any ejecta 
moving with a velocity of >1000 kms _1 should reach this 
over the timescales of the 2012b event. Any such ejecta will 
pass through the shocks, being decelerated and heated in 
the process. The fact that ejecta with velocities close to 10 4 
kms -1 is still observed in our last spectra indicates that ei- 
ther the photosphere is not at the interaction region - or that 
the cool dense shell/interaction region is clumpy, and that 
some of the fast ejecta can pass through it unhindered. The 
latter scenario could also help explain why even at the peak 
of the 2012b event there are weak broad absorptions. Such 
dumpiness is almost a certainty, given likely instabilities as- 
sociated with the the multiple shocks passing through the 
CSM due to past explosive events. 

An alternative arrangement to explain the simultaneous 
presence of high velocity absorption and interaction domi- 
nated line profil es could be the a disk o r torus-like geometry, 
as suggested by iLevesque et ail l|2012h . In such a scenario, 
most of the CSM lies in a disk, the polar axis of which is 
inclined relative to our line of sight. A fraction of the fast 
ejecta will impact the disk and give rise to the interaction, 
depending on the thickness of the disk and the precise emit- 
ting geometry. 

We suggest that the p ulsati onal-pair instability mod- 
els of IWooslev et all (|2007l ) and IChatzopoulos fc Wheeled 
l|2012l ) reproduce well the observed characteristics of SN 
2009ip, and that the high luminosity seen during the 2012b 
event is probably powered by colliding shells. The 2012a 
event could also be powered by interaction, alternatively it 
could be simply due to the internal energy released during 
an outburst. Either way, this prediction can be tested by 
observing the surviving progenitor star (or finding it has 
disappeared) once SN 2009ip emerges from behind the Sun 
in April 2013. Accurate measurements of the metallicity at 
the site of SN 2009ip would also be of interest, given the ex- 
pectation that pulsational-pair instability eruptions should 
only occur at low metallicity. 



7 CONCLUSIONS 

We have presented intensive UV, optical and NIR spec- 
troscopy and photometry of the interacting transient SN 
2009ip. The transient was intrinsically luminous at the peak 
of the 2012b event, but over the following weeks faded much 
more rapidly than is typically seen in CCSNe. From our 
photometry we have set a robust upper limit of <0.02 Moto 
the mass of any 56 Ni synthesized in SN 2009ip. This limit 
can be considered secure with respect to both 7-ray leakage 
and dust formation. Spectoscopically, the transient shows a 
blue continuum with narrow emission lines from H and He, 
before developing strong broad absorption and emission in 
both H and He, along with Na and Ca. At no point do we 
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see the nebular emission from synthesized material which is 
the hallmark of late stage evolution in CCSNe. 

From our modelling of the CSM, we find densities of 
10 9 -10 10 are sufficient to explain the narrow Balmer lines 
and the observed Balmer decrement. We also find that a 
velocity gradient is required in the CSM to reproduce the 
Ha to continuum ratio; such a gradient implies that the mass 
loss is eruptive rather than a steady wind. 

We see no evidence for core-collapse in our data set, and 
in particular do not observe the appearance of nebular emis- 
sion lines or detect a radioactive tail. We suggest that the 
properties of SN 2009ip can be explained with the efficient 
conversion of kinetic energy to optical luminosity via collid- 
ing shells, possibly resulting from a pair-instability driven 
outburst. 

A final question posed by SN 2009ip - if it is indeed not 
a core-collapse SN - is whether o ther interacting tran sients 
have been similarly misclassified. iDessart et al.l (|2009h sug- 
gested that SN 199 4W was not a CCSN, w hile the possibility 
was also raised by iKankare et~all l |2012h for SN 2009kn. If 
the progenitor of SN 2009ip is seen to survive this latest 
eruption, it will confirm that a high peak luminosity and 
the detection of high velocity ejecta cannot be regarded as 
definitive proof of core-collapse in interacting transients. 
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APPENDIX A: NLTE MODELLING FOR THE 
NARROW H LINES IN SN 2009IP 

We compute NLTE solutions for the ionization and exci- 
tation balance of hydrogen, taking all radiative processes 
due to the photospheric radiation field into account, along 
with all bound-bound collisional transitions. We use a 
model atom including fine structure with TV = 210 nl- 
levels (n = 1 — 20), with energy levels and A-values from 
NIST, bound-bound c ollision rates from I Johnson! (|l972f ) and 
Anderson et"afl i200Ch . photoionization cross sections from 

1 and and scaled to n ^ 2 ac- 



Verner et all dl996t) for 



cording to lRvbicki fc Lightman ( 1979), and radiative recom- 
bination rates from iBrocklehurstl (jl97l|). We do not include 
collisional ionization or its inverse process three-body re- 
combination, which start to become important only at high 
densities (n > 10 13 cm -3 ). 

We take radiative transfer effects into account by using 
escape probabilities. In the limit of small velocity gradients 
(dV/dr <C AVthermai/ Ar) , we use th e position-averaged es- 
cape probability in a uniform sphere (jOsterbrock fc Ferlandl 

mom 



Pirn) = — 



1 - 



2r 2 



8tt 3 / 2 Qi 



+ — 



9j 



2r 2 
AV tl 



Ar 



(Al) 



(A2) 



where i and j refer to the levels (i is the lower one here), 
gi and gj are the statistical weights, Aij is the Einstein co- 
efficient for spontaneous decay, Aij is the line wavelength, 
AVthermai is the thermal velocity dispersion 



AVthermai 



2fcT 



(A3) 



and Ar is the effective column length, given by 

Ar = 1 f Rwn ru{r) - nj (r)^dr (A4) 

ni{R in ) - nj(R in )f- J R . n g, 

where Ri on is the outer limit of the ionized region (computed 
iteratively) . 

In the limit of large velocity gradients (dV/dr 
AVthermai / Ar ) , we instead use the Sobolev approximation 
l|Sobolevlll960l : ICastorlll970h 



/?(7 



1 - e~ T 'j 



(A5) 



Since t^/t^ « AVthermai/ AVfiow, where AVfi ow ~ 
-Ar, rapid flows (AVfiow > AVthermai) have the effect 



9 Various flavors of escape probabilities may be used depend- 
ing on the assumed geometry and density profile of the emit- 
ting region as we ll as the details of the radiative trans- 
fer approximations jO stcrbrock & Fcrland 2006; Capriottil ll965l ; 
brake fc Ulrichl 1 1980 )7 but for our purposes trying to select the 
"best" one is not meaningful, as they only differ by a factor of a 
few at most. Note that we also use the escape probability at line 
center, rather than the line-profiled average expression, as these 
generally differ by less than a factor of two. 



of lowering the line optical depths, enhancing the escape 
probabilities. 

A priori, the velocity gradients in the CSM of SN2009ip 
are unknown. The gradient could quite possibly be impor- 
tant, since typical expansion velocities of the narrow lines in 
the spectra are of order AVfiow ~ 10 3 km s _1 , whereas the 
thermal velocities are of order AVthermai ~ 10 km s _1 (as- 
suming T~ 10 4 K). If the CSM was created by an eruption, 
the velocity field will be close to homology with a gradient 
of — t" 1 where t is the time since the eruption. If the 
CSM instead originates in a steady wind, and accelerations 
in the wind are small, the gradient can be insignificant over 
the emission region. 

We explore models with different assumptions about the 
velocity field. We introduce a default velocity gradient 
by taking the ratio of the typical CSM velocity (10 3 km s _1 ) 
to the system scale size (10 4 Rq), so 

^ =1.4.10- 7 s-\ (A7) 
dr o 

For an eruptive origin, this corresponds to t = 80 days. 
We later try models with 4^ = e^ n , for different values of 



NLTE solutions 

The bound-bound radiative rat es due to the external radi a- 
tion field is computed by (e.g.. ICastor]|l970l ; lElitzurlll982h 



Rij = BijJ Uij j3(nj) s 1 



(A8) 



where Bij is the Einstein coefficient for photoabsorption [i < 
j) and stimulated emission (i > j), J Vij is the mean intensity 
at the line frequency i/y (excluding any attenuation by the 
line), and P{nj) is the escape probability (described above). 
Close to the photosphere we have 



J„ = — / d<f> / B v sin 8 cos ( 
47r Jo Jo 

The photoionization rate for level i is 



Ti = 



I IV 



at.dv s 1 



(A9) 



(A10) 



where a l v is the ionization cross section for level i, v\ is the 
threshold frequency. The total photoionization rate is 



r = Vr ( i 

TIL 



UHI 



(AH) 



where nui = S»=i Ui ana * Ui ^ s * ne l eve l population of level 
i (computed iteratively). 

The ionization balance equation is, approximating the 
proton density to equal the electron density (n+ = n e ), 



Tn(l - x e ) = (a + a sti m) n x e 



(A12) 



where x e is the electron fraction (= n e /n) a is the effec- 
tive (including optical depth effects) spontaneous radiative 
recombination coefficient 



a = ^a i( S(ri) 



(A13) 



where P(ji) is computed with the static escape probability 
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formalism (Eg. IA1[) . and a s 
tion coefficient 



is the stimulated recombina- 



atsti 



1=1 v 



4nJ » i -(h„-hv>„)/kT 



hv 



dv 



(A14) 

Eq. IA12I is solved by the Newton- Raphson method. 

For a given solution to x e (and thereby for n+ = n e = 
nx e and nni — n (1 — x e )), we then solve the set of statis- 
tical equilibrium equations for the excited levels of HI. For 
level i: 

[ai^Ti) + a£ tim ] "a + ^' l C 3ine + Ajifiji + Rji] (A15) 



r, + 



(A16) 



where Qj is the specific radiative recombination coefficient, 
dj are the bound-bound collision coefficients, and the other 
symbols are defined above. Since the escape probabilities fiij 
depend on the level populations themselves, the system is 
non-linear and is solved iteratively. 

Together with the number conservation equation (we ig- 
nore contributions by He and metals to n e , since this contri- 
bution is small unless very high ionization states are reached 
for the metals) 



iihi + n e = n 



(A17) 



we then iteratively solve for ionization balance, level popu- 
lations, and line optical depths, for given values of n and T„ 
(T is taken as = T,). This allows us to compute line optical 
depths and emissivity ratios for models in the Sobolev limit, 
where these quantities only depend on local conditions 

If we assume a density profile for the CSM, we can also 
compute solutions for steady flows (dV/dr = 0), as well 
as integrated quantities such as total line luminosities and 
continuum optical depths in all the models. The size of the 
ionized region can be estimated by equating the number of 
absorbed continuum photons to the number of recombina- 
tions 



-S>j(flio„) 



dv 



4-7rr n(r) adr 



(A18) 

which can be solved for Ri on given a density profile n(r). 
The continuum optical depths are computed by Eq. IA21I 

For small enough densities, there is no solution as re- 
combinations are too slow to bound an H n region, and the 
ionized region in then formally infinite (in actuality it is 
bounded by time travel limits, or when it encounters the 
interstellar medium.) 

Given the solution to Rion, the electron scattering op- 
tical depth is 



,(r)dr 



which for n e (r) = n e (R/R in ) p (p 2) equals 

,p-i 



M 

p-i 



(A19) 



(A20) 



where a e3 = 6.65 • 10 



g . For the level populations, 



it is a reasonable approximation that flj(r) oc r p , since re- 
combination rates scale with n 2 . ~ n 2 (through the dilution 
of the photospheric radiation field for photoexcitation and 
through the n 2 dependence for recombination. Thus, for the 
photoionization optical depths 



Oi{vl)ni(R in )Ri, 



(fc) 



2p-l 



2p 



(A21) 



To clarify, we thus compute the NLTE solutions ex- 
plicitly for the density and radiation field specified at the 
inner region, and assume scalings of the solutions at larger 
radii according to n e oc r~ p and rn oc r~ 2p . In all the solu- 
tions presented in the paper, we use p=2, corresponding to 
a steady mass loss. 
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Table Bl: Optical photometry of SN 2009ip in the Landolt system, un- 
certainties are given in parentheses. 



Date 


MJD 




U 




B 






V 






R 








I 




Instrument 


09/10/12 


56209 


04 


13.00 (0.04) 


13 


87 


(0 


02) 


13 


90 


(0 


01) 


13 


75 


(0 


01) 


13 


68 


(0 


01) 


EFOSC2 


10/10/12 


56210 


27 


13.04 (0.02) 


13 


93 


(0 


01) 


13 


93 


(0 


01) 


13 


78 


(0 


01) 


13 


69 


(0 


01) 


EFOSC2 


15/10/12 


56215 


24 


13.15 (0.02) 


14 


01 


(o 


01) 


13 


97 


(0 


01) 


13 


79 


(0 


01) 


13 


72 


(0 


01) 


EFOSC2 


16/10/12 


56216 


02 


13.22 (0.01) 


14 


09 


(o 


01) 


14 


06 


(o 


01) 


13 


87 


(0 


01) 


13 


79 


(0 


01) 


EFOSC2 


16/10/12 


56216 


71 


13.34 (0.05) 


14 


16 


(o 


05) 


14 


03 


(o 


05) 


















UVOT 


16/10/12 


56216 


77 


- 


14 


28 


(o 


05) 




- 


















- 




UVOT 


18/10/12 


56218 


01 


- 


14 


40 


(o 


02) 


14 


17 


(o 


02) 


13 


97 


(0 


02) 


13 


92 


(0 


04) 


PROMPT 


19/10/12 


56219 


02 


- 




- 






14 


21 


(0 


01) 


14 


04 


(0 


03) 


13 


95 


(0 


03) 


PROMPT 


18/10/12 


56219 


32 


13.82 (0.05) 


14 


37 


(o 


05) 


14 


20 


(0 


05) 














- 




UVOT 


20/10/12 


56220 


89 


13.84 (0.03) 


14 


45 


(o 


06) 


14 


28 


(o 


04) 


14 


09 


(0 


04) 


13 


99 


(0 


05) 


RATCAM 


20/10/12 


56220 


96 


13.88 (0.04) 




- 






























UVOT 


20/10/12 


56221 


04 


13.90 (0.05) 


14 


47 


(o 


05) 


14 


30 


(o 


05) 


















UVOT 


20/10/12 


56221 


18 


13.80 (0.04) 


14 


51 


(o 


01) 


14 


34 


(o 


01) 


14 


15 


(0 


01) 


14 


02 


(0 


01) 


EFOSC2 


21/10/12 


56221 


91 


13.89 (0.02) 


14 


51 


(o 


02) 


14 


35 


(o 


01) 


14 


15 


(0 


01) 










RATCAM 


22/10/12 


56222 


78 


13.93 (0.05) 


14 


61 


(o 


05) 


14 


45 


(o 


05) 


















UVOT 


22/10/12 


56222 


92 


14.07 (0.10) 


14 


68 


(o 


02) 


14 


41 


(o 


02) 


14 


19 


(0 


02) 


14 


15 


(0 


02) 


RATCAM 


23/10/12 


56223 


07 


- 




- 






14 


40 


(o 


01) 


14 


22 


(0 


01) 


14 


15 


(0 


04) 


PROMPT 


23/10/12 


56223 


25 


14.13 (0.03) 


14 


59 


(o 


02) 


14 


44 


(o 


01) 


14 


24 


(0 


01) 


14 


12 


(0 


02) 


EFOSC2 


24/10/12 


56224 


87 


14.20 (0.02) 


14 


74 


(o 


01) 


14 


48 


(o 


02) 


14 


29 


(0 


01) 


14 


22 


(0 


05) 


RATCAM 


24/10/12 


56224 


99 


14.26 (0.06) 


14 


66 


(o 


03) 


14 


49 


(o 


06) 


















UVOT 


25/10/12 


56225 


02 


- 


14 


79 


(o 


02) 


14 


54 


(o 


05) 


14 


29 


(0 


03) 


14 


22 


(0 


02) 


PROMPT 


24/10/12 


56225 


12 


14.27 (0.05) 




- 








- 


















- 




UVOT 


25/10/12 


56225 


88 


14.33 (0.09) 


14 


86 


(o 


02) 


14 


65 


(o 


01) 


14 


41 


(0 


01) 


14 


34 


(0 


03) 


RATCAM 


26/10/12 


56226 


73 


14.62 (0.06) 


14 


92 


(o 


05) 


14 


69 


(o 


06) 


















UVOT 


28/10/12 


56228 


86 


14.97 (0.06) 




- 






























UVOT 


29/10/12 


56229 


05 


- 


15 


42 


(o 


05) 


15 


02 


(o 


03) 


14 


74 


(0 


04) 


14 


61 


(0 


04) 


PROMPT 


28/10/12 


56229 


21 


15.03 (0.06) 


15 


23 


(o 


06) 


15 


00 


(o 


07) 


















UVOT 


30/10/12 


56230 


87 


15.38 (0.06) 


15 


57 


(o 


05) 


15 


13 


(o 


06) 


















UVOT 


31/10/12 


56231 


03 


- 


15 


49 


(o 


03) 




- 






















PROMPT 


01/11/12 


56233 


35 


15.25 (0.06) 


15 


39 


(o 


06) 


14 


97 


(o 


06) 


















UVOT 


03/11/12 


56234 


98 


15.17 (0.05) 


15 


22 


(o 


05) 


14 


78 


(o 


06) 


















UVOT 


05/11/12 


56236 


09 


14.87 (0.01) 


15 


21 


(o 


01) 


14 


92 


(o 


06) 


14 


71 


(0 


01) 


14 


36 


(0 


04) 


EFOSC2 


05/11/12 


56236 


71 


15.02 (0.05) 


15 


17 


(o 


06) 


14 


81 


(o 


05) 


















UVOT 


06/11/12 


56237 


07 


- 


15 


32 


(o 


03) 


























PROMPT 


07/11/12 


56238 


03 


- 


15 


37 


(o 


04) 


15 


07 


(o 


03) 


14 


84 


(0 


03) 


14 


68 


(0 


03) 


PROMPT 


07/11/12 


56238 


55 


15.39 (0.05) 


15 


40 


(o 


04) 


15 


15 


(o 


07) 


















UVOT 


08/11/12 


56239 


02 


- 


15 


56 


(o 


02) 




- 


















- 




PROMPT 


08/11/12 


56239 


14 


15.33 (0.01) 


15 


56 


(o 


01) 


15 


26 


(o 


01) 


14 


96 


(0 


02) 


14 


80 


(0 


02) 


EFOSC2 


09/11/12 


56240 


55 


15.68 (0.05) 


15 


74 


(o 


05) 


15 


29 


(o 


06) 


















UVOT 


11/11/12 


56242 


04 


- 


15 


99 


(o 


05) 


15 


50 


(o 


02) 


15 


21 


(0 


02) 


15 


03 


(0 


04) 


PROMPT 


11/11/12 


56242 


79 


16.11 (0.07) 


16 


14 


(o 


06) 


15 


57 


(o 


06) 


















UVOT 


11/11/12 


56242 


80 


16.08 (0.05) 


16 


09 


(o 


02) 


15 


60 


(o 


05) 


15 


26 


(0 


07) 


15 


03 


(0 


05) 


RATCAM 


12/11/12 


56243 


04 


- 


16 


13 


(o 


04) 


15 


59 


(o 


02) 


15 


24 


(0 


03) 


15 


11 


(0 


04) 


PROMPT 


13/11/12 


56244 


04 


- 


16 


17 


(o 


02) 


15 


67 


(o 


02) 


15 


35 


(0 


02) 


15 


15 


(0 


04) 


PROMPT 


13/11/12 


56244 


92 


16.39 (0.07) 


16 


24 


(o 


06) 


15 


71 


(o 


06) 














- 




UVOT 


13/11/12 


56244 


82 


16.281 (0.09) 


16 


30 


(o 


06) 


15 


69 


(o 


02) 


15 


35 


(0 


08) 


15 


22 


(0 


03) 


RATCAM 


14/11/12 


56246 


11 


- 


















15 


37 


(0 


02) 






- 




EFOSC2 


15/11/12 


56246 


17 




16 


29 


(o 


03) 


15 


70 


(o 


07) 


15 


42 


(0 


03) 


15 


22 


(0 


05) 


PROMPT 


15/11/12 


56247 


41 


16.52 (0.08) 


16 


41 


(o 


05) 


15 


81 


(o 


07) 


















UVOT 


17/11/12 


56248 


54 


16.58 (0.06) 


16 


43 


(o 


06) 


15 


85 


(0 


07) 


















UVOT 


18/11/12 


56249 


08 




16 


43 


(o 


03) 


15 


80 


(o 


03) 


15 


48 


(0 


02) 


15 


27 


(0 


04) 


PROMPT 


19/11/12 


56250 


54 


16.87 (0.09) 


16 


59 


(o 


07) 


15 


96 


(o 


07) 


















UVOT 


21/11/12 


56252 


13 




16 


76 


(o 


03) 


16 


04 


(o 


06) 


15 


69 


(0 


04) 


15 


45 


(0 


06) 


PROMPT 


21/11/12 


56252 


16 


17.12 (0.02) 


16 


77 


(o 


01) 


16 


21 


(o 


06) 


15 


73 


(0 


04) 


15 


48 


(0 


02) 


EFOSC2 


21/11/12 


56253 


08 


17.15 (0.09) 


16 


82 


(o 


08) 


16 


15 


(o 


08) 


















UVOT 


23/11/12 


56254 


58 


17.51 (0.09) 


16 


87 


(o 


07) 


16 


29 


(o 


07) 


















UVOT 


25/11/12 


56256 


14 




16 


95 


(o 


05) 


16 


25 


(o 


04) 


15 


81 


(0 


04) 


15 


52 


(0 


05) 


PROMPT 


25/11/12 


56257 


17 


17.64 (0.08) 


17 


04 


(o 


07) 


16 


38 


(o 


09) 


















UVOT 
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27/11/12 


56259.03 


17.66 (0.08) 


17.19 


(0.08) 


16.47 


(0.07) 


- 


- 


UVOT 


29/11/12 


56260.83 


17.64 (0.09) 


17.19 


(0.07) 


16.52 


(0.07) 


- 


- 


UVOT 


30/11/12 


56261.15 


- 


17.06 


(0.04) 


16.37 


(0.03) 


- 


- 


PROMPT 


01/12/12 


56262.70 


17.62 (0.08) 


17.18 


(0.07) 


16.47 


(0.07) 


- 


- 


UVOT 


03/12/12 


56264.08 


- 


17.22 


(0.04) 


16.40 


(0.01) 


15.93 (0.02) 


15.62 (0.05) 


PROMPT 


03/12/12 


56264.82 


17.46 (0.10) 


17.14 


(0.07) 


16.48 


(0.04) 


15.86 (0.08) 


15.70 (0.05) 


RATCAM 


03/12/12 


56265.06 


17.66 (0.10) 


17.31 


(0.06) 


16.54 


(0.04) 


16.06 (0.03) 


15.66 (0.05) 


EFOSC2 


03/12/12 


56265.11 


17.82 (0.09) 


17.37 


(0.08) 


16.69 


(0.07) 


- 


- 


UVOT 


05/12/12 


56266.65 


18.01 (0.12) 


17.68 


(0.09) 


16.72 


(0.08) 


- 


- 


UVOT 


07/12/12 


56269.16 


18.59 (0.17) 


18.16 


(0.11) 


17.04 


(0.11) 


- 


- 


UVOT 


09/12/12 


56270.85 


18.87 (0.24) 


18.31 


(0.14) 


17.32 


(0.13) 


- 


- 


UVOT 


11/12/12 


56272.08 


- 


18.34 


(0.04) 


17.30 


(0.04) 


16.62 (0.03) 


16.32 (0.09) 


PROMPT 


11/12/12 


56272.76 


18.97 (0.21) 


18.52 


(0.14) 


17.46 


(0.14) 






UVOT 


11/12/12 


56273.05 












16.72 (0.05) 




EFOSC2 


13/12/12 


56274.62 


19.06 (0.42) 


18.86 


(0.33) 


18.06 


(0.30) 






UVOT 


14/12/12 


56275.07 




19.03 


(0.06) 


18.00 


(0.05) 


17.22 (0.04) 


16.94 (0.07) 


PROMPT 


15/12/12 


56276.73 




19.30 


(0.22) 


18.53 


(0.21) 






UVOT 


17/12/12 


56278.05 




19.16 


(0.10) 


18.35 


(0.14) 


17.46 (0.11) 


17.24 (0.11) 


PROMPT 


20/12/12 


56282.09 


19.452 (0.10) 


19.22 


(0.09) 


18.27 


(0.06) 


17.42 (0.05) 


17.15 (0.10) 


EFOSC2 


29/12/12 


56290.07 




>18.0 


18.55 


(0.11) 






PROMPT 
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Table B2. UV photometry of SN 2009ip from Swift+WOT, together with errors in parentheses. All values are in the UVOT system. 



Date MJD uvw2 uvm2 uvwl 



lb/ 1U/12 


56216 


71 


13 


39 


(0 


05) 


13 


12 


(0 


05) 


13 


33 


(° 


05) 


i z" 1 / "i r\ /"io 

16/10/12 


56216 


87 


13 


54 


(0 


05) 


















16/10/12 


56216 


90 










13 


29 


(0 


06) 










18/10/12 


56219 


32 


14 


20 


(0 


05) 


13 


97 


(0 


06) 


13 


79 


(0 


06) 


20/10/12 


56220 


93 


14 


58 


(0 


05) 


















20/10/12 


56221 


04 


14 


59 


(0 


06) 


14 


39 


(0 


06) 


14 


10 


(0 


06) 


22/10/12 


56222 


78 


14 


90 


(0 


06) 


14 


55 


(0 


06) 


14 


26 


(0 


06) 


22/10/12 


56222 


98 










14 


56 


(0 


07) 










no /-in / 1 o 

22/ 10/12 


56222 


99 


14 


70 


(0 


05) 


















n yi / 1 n / 1 o 

24/ 10/12 


56224 


99 


15 


29 


(0 


07) 


14 


93 


(0 


06) 


14 


48 


(0 


06) 


24/10/12 


56225 


00 










15 


15 


(0 


07) 










24/10/12 


56225 


14 


15 


31 


(0 


06) 


















26/10/12 


56226 


73 


15 


54 


(0 


07) 


15 


26 


(0 


07) 


14 


85 


(0 


06) 


28/10/12 


56229 


21 


15 


98 


(0 


08) 


15 


61 


(0 


07) 


15 


39 


(0 


07) 


30/10/12 


56230 


87 


16 


51 


(0 


09) 


16 


00 


(0 


05) 


15 


59 


(0 


07) 


01/11/12 


56233 


35 


16 


64 


(0 


08) 


16 


19 


(0 


07) 


15 


76 


(0 


07) 


02/11/12 


56234 


28 


16 


64 


(0 


08) 


16 


27 


(0 


04) 










03/11/12 


56234 


98 


16 


68 


(0 


08) 


16 


20 


(0 


06) 


15 


64 


(0 


06) 


05/11/12 


56236 


71 


16 


68 


(0 


07) 


16 


42 


(0 


06) 


15 


65 


(0 


06) 


/""IT / 1 1 / 1 O 

07/ 11/12 


56238 


55 


16 


51 


(0 


09) 


16 


20 


(0 


07) 


15 


58 




08) 


09/11/12 


56240 


55 


17 


19 


(0 


06) 


16 


74 


(0 


06) 


16 


06 


(0 


08) 


11/11/12 


56242 


79 


17 


57 


(0 


10) 


16 


91 


(0 


09) 


16 


45 


(0 


08) 


13/11/12 


56244 


92 


17 


81 


(0 


10) 


17 


51 


(0 


10) 


17 


08 


(0 


16) 


15/11/12 


56247 


41 


18 


00 


(0 


15) 


17 


70 


(o 


13) 


17 


26 


(0 


10) 


17/11/12 


56248 


54 


18 


15 


(0 


12) 


17 


80 


(0 


11) 


17 


28 


(0 


10) 


19/11/12 


56250 


54 


18 


43 


(0 


11) 


17 


98 


(0 


12) 


17 


35 


(0 


12) 


21/11/12 


56253 


08 


18 


94 


(0 


17) 


18 


68 


(0 


15) 


17 


57 


(0 


11) 


23/11/12 


56254 


58 


19 


19 


(0 


16) 


18 


75 


(0 


10) 


17 


75 


(0 


11) 


25/11/12 


56257 


17 


19 


45 


(0 


22) 


19 


06 


(0 


17) 


18 


07 


(0 


10) 


27/11/12 


56259 


03 


19 


68 


(0 


22) 


19 


26 


(0 


18) 


18 


22 


(0 


11) 


29/11/12 


56260 


83 


20 


08 


(0 


26) 


19 


14 


(0 


13) 


18 


36 


(0 


11) 


01/12/12 


56262 


70 


20 


10 


(0 


26) 


19 


59 


(0 


17) 


18 


72 


(0 


13) 


03/12/12 


56265 


11 


20 


10 


(0 


26) 


19 


56 


(0 


20) 


18 


68 


(0 


14) 



Table B3. NIR photometry of SN 2009ip in the 2MASS system. Errors are given in parentheses. 



Date 


MJD 








J 






H 




K 


Instrument 


12/10/12 


56212 


24 


13 


98 


(0 


07) 


13 


78 


(0 


04) 






REMIR 


19/10/12 


56330 


16 










13 


49 


(0 


07) 






LUCIFER 


20/10/12 


56220 


83 


13 


99 


(0 


08) 


13 


67 


(0 


04) 


13.27 


(0.08) 


NOTCAM 


21/10/12 


56221 


02 


14 


19 


(0 


13) 














REMIR 


22/10/12 


56223 


19 


14 


17 


(0 


02) 


13 


71 


(0 


04) 


13.66 


(0.13) 


SOFI 


26/10/12 


56226 


06 


14 


40 


(0 


12) 














REMIR 


05/11/12 


56237 


05 


14 


41 


(0 


08) 


14 


22 


(0 


11) 


13.94 


(0.12) 


SOFI 


11/11/12 


56242 


18 


15 


18 


(0 


08) 


Id 


97 


(0 


10) 






REMIR 


13/11/12 


56245 


06 


14 


95 


(0 


07) 


14 


75 


(0 


08) 


14.57 


(0.06) 


SOFI 


04/12/12 


56266 


02 


15 


33 


(0 


04) 


15 


09 


(0 


08) 


14.91 


(0.05) 


SOFI 


21/12/12 


56283 


08 


16 


88 


(0 


11) 


16 


62 


(0 


08) 


16.29 


(0.09) 


SOFI 



© 0000 RAS, MNRAS 000, 000-000 



